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Abstract. Neutrinos are fundamental particles ubiquitous in the Universe and whose 
properties remain elusive despite more than 50 years of intense research activity. This 
review illustrates the importance of solar neutrinos in Astrophysics, Nuclear Physics, 
and Particle Physics. After a description of the historical context, we remind the reader 
of the noticeable properties of these particles and of the stakes of the solar neutrino 
puzzle. The Standard Solar Model triggered persistent efforts in fundamental Physics 
to predict the solar neutrino fluxes, and its constantly evolving predictions have been 
regularly compared to the detected neutrino signals. Anticipating that this standard 
model could not reproduce the internal solar dynamics, a Seismic Solar Model was 
developed which enriched theoretical neutrino flux predictions with in situ observation 
of acoustic and gravity waves propagating in the Sun. This seismic model contributed 
to the stabilization of the neutrino flux predictions. This review reminds the main 
historical steps, from the pioneering Homestake mine experiment and the GALLEX- 
SAGE experiments capturing the first pp neutrinos. It emphasizes the importance 
of the Superkamiokande and SNO detectors. Both experiments demonstrated that 
the solar-emitted electronic neutrinos are partially transformed into other neutrino 
flavors before reaching the Earth. This sustained experimental effort opens the door 
to Neutrino Astronomy, with long-base lines and underground detectors. The success 
of BOREXINO in detecting the ^ Be neutrino signal alone instills confidence in the 
physicists ability to detect each neutrino source separately. It justifies the building of a 
new generation of detectors to measure the entire solar neutrino spectrum with greater 
detail, as well as supernova neutrinos. A coherent picture emerged from neutrino 
physics and helioseismology. Today, new paradigms take shape in these two fields: the 
neutrinos are massive particles, but their masses are still unknown, and the research on 
the solar interior is focusing on the dynamical aspects and on signature of dark matter. 
The magnetic moment of the neutrino begins to be an actor of stellar evolution. The 
third part of the review is dedicated to this prospect. The understanding of the crucial 
role of both rotation and magnetism in solar physics benefit from SoHO, SDO, and 
PICARD space observations, and from new prototype like GOLF-NG. The magneto- 
hydrodynamical view of the solar interior is a new way of understanding the impact 
of the Sun on the Earth environment and climate. For now, the particle and stellar 
challenges seem decoupled, but this is only a superficial appearance. The development 
of asteroseismology — with the COROT and KEPLER spacecrafts — and of neutrino 
physics will both contribute to improvements in our understanding of, for instance, 
supernova explosions. This shows the far-reaching impact of Neutrino and Stellar 
Astronomy. 

PACS numbers: solar neutrinos, neutrino properties, global helioseismology, 
local helioseismology, internal solar rotation, internal solar magnetic fields, Sun- 
Earth relationship. 
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1. Introduction 

The Sun is an outstanding and permanent source of neutrinos of different energy. 
Understanding this source is crucial to different scientific disciphnes, and a parallel 
progress on neutrino properties is necessary to allow important developments in Neutrino 
Astronomy. The Sun is the nearest neutrino source in the Universe (see table 1 
and Figure 1) and the first source in term of detectable neutrinos (it produces about 
6.7 X 10^*^z^/cm^/s), hence its importance. The strongest neutrino source, the cosmic 
neutrino background radiation, has a flux of 10^^z//cm^/s, but the energy of these 
neutrinos is too low (about 10"'^ — 10~^ eV) to currently allow a direct detection. 

In the late sixties, it became clear that we did not understand the solar-neutrino 
properties. From there, it took about half a century to provide solutions to the 
"solar neutrino puzzle", or deficit of neutrinos detected on Earth compared with the 
"theoretical estimate" of the emitted fiuxes. We now reached a coherent picture of 
this problem. We know unambiguously that neutrinos have a mass, unlike what was 
thought for a long time. However a lot of questions remain unanswered regarding the 
real nature of the neutrinos, their magnetic interaction with plasma, with the internal 
solar dynamics, and their connection to the solar activity. 

During the last twenty years, helioseismology provided Astrophysics with an 
opportunity to describe with an unprecedented accuracy the different solar neutrino 
sources and to look for coherence between the two existing probes of the deep solar 
core. In parallel, this discipline stimulated the development of research on the internal 
Stellar Dynamics. Today an impressive agreement between "helioseismic" neutrino fiux 
predictions and neutrino detections by the existing detectors has been reached. Thanks 
to this agreement, it becomes possible to describe the solar neutrino sources with an 
accuracy equivalent to the one reached from particle accelerators or nuclear reactors 
producing neutrinos or antineutrinos. Combined with the improvements in neutrino 
detectors, these facts favor the development of Neutrino Astronomy. Considering the 
vitality of the two disciplines (including asteroseismology, the seismology of stars), 
important discoveries are bound to be made in the upcoming decades. 

With hindsight, it is clear that using the second best source of neutrinos in 
the Universe, the Sun, was especially judicious because it allowed great advances in 
numerous fields of physics which contribute to a better description of both stars and 
neutrino properties. This review details the main aspects of this research (sections 3 to 
7). Section 8 details currently unanswered and interesting questions justifying upcoming 
developments, including the question of dark matter. First, we start with an historical 
review of the genesis of this research field (section 2). 

2. Historical context: the genesis of this field 

Pauh was still young, on December 4*^^, 1930, when he suggested the existence of a 
"small neutral particle" to interpret the continuum energy spectrum of the electrons 
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pp chain: 98.8% of the total energy produced by the present Sun 

p + p— 7'D + e"'" + fe called pp neutrinos 
(0.25%) :p + p + e~— )-D + i/e called pep neutrinos 
pp I chain (86%) : ^He +^ He -^"^ He + 2p Qes = 26.2 MeV 
pp n chain : (14%) : ^He +4 He ^'^ Be + 7 
^Be + e~ Li + called ^ Be neutrinos 
^Li + p ^ 2'^He Qeff = 25.66 McV 
pp III chain: (0.2%): '^Bc + p B Qcs = 19.17McV 
®B Be* +e+ + Uc ^Be* 2'^Hc called neutrinos 
pp IV chain : (0.002%) p +3 He i/e + e+ He (hep neutrinos) 

CNO cycle: 1.2% of the total energy produced by the present Sun 

CNO I cycle: ^^N -^^^ C + e+ + called ^^N neutrinos 

150 -^^^ N + e+ + i/e called ^^O neutrinos 
CNO II cycle : ^^F ^^'^ O + e+ + z/e called ^'^F neutrinos 
Qeff CNO = 26.73 MeV 




Figure 1. Energy dependence of the neutrino sources and the range of energy covered 
by the neutrino detectors. 
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in the /3 radioactivity, through a letter pubhcly read at the Tubingen conference. It 
is difficuh to beheve that the existence of this "neutrino" particle was proposed before 
the discovery of the neutron (Chadwick 1932). This brilliant idea was immediately 
supported by Fermi in 1931 as an important piece in the weak interaction theory, and 
he coined the name of "the little neutral" . Nevertheless, the first demonstration of its 
existence by F. Reines and C. Cowan, near the reactor of Savannah River, came only 
25 years later, in 1956. 

Pontecorvo, in 1946, suggested to use the Sun as an important source of neutrinos, 
and Crane (1948) showed the potentiality of neutrino physics. In parallel, the knowledge 



on stellar evolution was steadfastly progressing. Eddington (1926) had already 
understood that the continuous brightness of the stars was of nuclear origin. Then, 
Gamow, Weizsacker, and Bethe solved the conundrum of these nuclear interactions: the 
competition between the Maxwellian tail distribution and the presence of the Coulomb 
barrier, and the existence of different chains of reactions (see table 1 and section 4). 
The number of neutrinos produced in the Sun can be directly deduced from the fact 
that 4 protons convert hydrogen to helium while releasing 2 neutrinos and some energy. 
The neutrino flux is obtained from the ratio of the solar luminosity to the quantity of 
energy released. 

4p He + 2e+ + 2z/e + 26.20MeV - 

where 26.20 MeV is the most common effective energy and Ey has a mean value of 
0.26 MeV (in 85% of the cases, see tables 1 and 2). Based on this reaction, about 7 
X 10^'^neutrinos/ cm^/s escape from the Sun. Undoubtedly these developments fascinated 
the young R. Davis, who received a formation in physical chemistry and was looking 
for a good research topic, when he joined the Brookhaven National Laboratory in the 
aftermath of the second world war. 

After several years spent on looking for the best detection technique, Davis decided 
in 1951 to follow the idea of Pontecorvo and to try capturing neutrinos through the 
reaction: '^''Cl + Vc -^^'^ Ar + e~, with a decay back to ^^Cl by the inverse capture 
process with a half-life of 35 days. There were at that time two possible sources of 
detectable neutrinos: the nuclear fission reactors and the Sun. Therefore Davis first 
tried to validate this method with a 3900-liter tank of carbon-tetrachloride (CCI4) as the 
target material, located near the Brookhaven nuclear reactor. Back then, the fact that 
only antineutrinos, and no neutrinos, were produced was not deemed a good enough 
reason to pass this opportunity. The reactor was considered not powerful enough to 
produce any neutrino signal. In parallel, the background noise resulting from cosmic 
rays was clearly identified, and a first experimental upper limit for the solar neutrino 
flux was set at 40000 SNU ( |Davis| |1955D (where 1 SNU, the Solar Neutrino Unit, is 
defined as 10"'^® captures per target atom per second). Even at the time, this result 
appeared uninformative. This first attempt also highlighted the necessity to recourse to 
underground experiments. Davis decided to repeat this experiment twice at Savannah 
River, the most powerful site of nuclear reactors in the world, first with the same quantity 
of carbon-tetrachloride, and then with a 11400-liter tank. The absence of any reactor 
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neutrino signal lowered the upper limit for neutrino capture to a factor 20 below the 



antineutrino capture rate. The first detection of a free antineutrino by Cowan et al. 



(1956); Reines et al. (1960), using a different technique, coupled with the results of 
Davis, showed that the neutrino was not its own antiparticle (see Davis (2002)). They 
will win the Nobel prize in 1995 for this discovery. See section 3 for a more general 
consideration on neutrino properties). 

The Sun produces a great many neutrinos. Unfortunately these neutrinos come 
mainly from the fundamental transformation of a proton into a neutron: p + p — ?• 
2D + i/g + e"*" (where D is a deuteron), with a resulting neutrino energy lower than 0.42 
MeV (see table 1 and figure 1). With such low energy these neutrinos were not accessible 
to the detector of R. Davis, since the energy threshold of his experiment was 0.814 
MeV. At that time, the ppl chain converting hydrogen into helium through (^He,^ He) 
interaction was considered as the main actor of the proton-neutron transformation. But 
when Holmgren and Johnston reported in 1958 that the reaction rate of (^He+^He) was 
1000 times greater than previously thought, the production of ''Be was then increased 
by a factor 7000, immediately raising the odds of a ^Be-proton interaction producing 
'^B with emission of an energetic neutrino, and also the odds of an electronic capture 
producing ^Li and another neutrino. Both reactions are energetic sources of neutrinos 
and the resulting prospect of several neutrino captures per day was a strong motivation 
for building a 3800-liter tank of perchloroethylene (C2CI4) in the Barberton Limestone 
mine in Ohio. This solution of perchloroethylene was widely used in dry cleaner shops 
and less dangerous for the health. Unfortunately, once more, Davis did not detect any 
neutrino signal. 

The different reaction cross sections were measured. Kavanagh (1960) measured 
the ^Be(p, 7) reaction rate and it appeared that this interaction cross-section was much 
smaller than the electron capture cross-section. Due to the energy range of the neutrinos 
produced by the two reactions (0 to 14.1 MeV for neutrinos associated with ^B and 
two narrow lines of respectively 0.38 and 0.86 MeV for those associated with ^Li), the 
likelihood of detecting solar neutrinos dropped, as written by Reines: "The probability 
of a negative result even with detectors of thousands or possibly hundreds of thousands 
of gallons of perchloroethylene tends to dissuade experimentalists from making the 
attempts" . 

However, Bahcall (1962) calculated the electron capture rate and the capture rate 
of *B and showed that the latter rate was 20 times larger than previously thought. 
Therefore there was a real hope of detecting solar neutrinos: the theory became a 
guidance for determining the first neutrino fluxes (Bahcall et al. , 1963) and the size of 
the detector (Bahcall, 1964a) : for instance, 378000 liters of perchloroethylene could 
detect 4 to 11 ^^Ar atoms per day (about 28 SNU). In parallel, the work on muon 
cosmic-ray at Barberton was useful to determine the depth of the set up needed to kill 
the main source of noise (Davis, 1964) Davis, in collaboration with Bahcall, put a tank 
in the Homestake Gold Mine in Dakota at the -1478 m level of the mine. When the 
measurements started in 1968, it was clear that the main contributor to the detected 
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neutrinos was the proton capture rate on ^Be. The neutrino flux predictions (Bahcall 



& Shaviv, 1968) included an error bar, due to the uncertainties on the ingredients of 



the solar model, of about 60%. In the following decades, as shown in this review, John 
Bahcall turned out to be the most active theoretician sustaining and stimulating work 
on new neutrino detection techniques, thanks to his extremely well documented neutrino 
flux predictions. He is at the root of the solar neutrino puzzle. 

In 1968, R. Davis and collaborators derived an upper limit of 3 SNU for the solar 
neutrino emission (Davis, Harmer & Hoffman, 1968) To obtain this result, he had first to 
use some ^^Ar in order to demonstrate his ability to extract argon atom at the 95% level. 
In 1971, he announced a detected flux of 1.5 ± 1 SNU. This value has not significantly 
changed over time, and except a reduction in the error bar partly due to the integration 
over a longer timespan, the statistical accuracy of each run remained very low (about 7 
events per month on average). In the meantime, J. Bahcall, N. Bahcall, & G. Shaviv, in 
a letter to Physical Review Letters in 1968, produced different models of the Sun and 
favored a model leading to a chlorine neutrino flux prediction of 7.7 SNU, much more 
than the detected flux. Immediately, Pontecorvo mentioned the possibility of oscillation 
between neutrino flavors (Pontecorvo, 1968). 1968 was undeniably a very inspiring year 
for neutrino physics. For his lifetime efforts, R. Davis received the Nobel prize in 2002, 
34 years later. This Nobel prize was shared with M. Koshiba for the outstanding success 
of the Kamiokande and Superkamiokande collaborations (discussed later in this review). 
The detection of atmospheric and solar neutrinos contributed to the evidence of neutrino 
oscillation (see below). 

The prediction of solar neutrino fluxes from a solar model results from complex 
calculations, and J. Bahcall played a major role in this field. Such predictions require 
a good estimate of the nuclear reaction rates listed in tables 1 or 2, a good estimate 
of the neutrino cross sections (Bahcall, 1964b, 1989), a good knowledge of the Sun, 
and a strong understanding of the neutrino properties. Nevertheless, convergence of 
the various estimates of the emitted neutrino flux was rather quickly achieved, and 
the predictions have not changed by more than a factor of 2 (for neutrinos) once 
the first period of "instability" settled (see all the publications dedicated to this topic 
between 1985 and 1998, and mentioned in the present article). However, the deficit of 
neutrino detection compared to the prediction, coined the neutrino puzzle as early as the 
seventies, clearly meant that a great deal of interesting physics was yet to be discovered: 
new properties of the neutrinos and/or new insight in the solar (stellar) physics beyond 
the standard framework. Since this standard framework (ignoring, among others, the 
role of the magnetic field) could not reproduce the observations, this raised the question 
of whether or not the Sun should be treated as a magnetic star. That means including 
a lot of processes that are potentially present in the Sun, such as mixing in the solar 



core, rotation, magnetic field ( 


Shaviv & Salpeter 


1971 


Schatzman 


al. 


1981 


Roxburgh , 


1985 


Gough & Thompson, 


1990 


) but that are i 



without constraints in the radiative zone. This period was very stimulating for all the 
research fields linked to the neutrino puzzle. Variations by up to 30% in the prediction 
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of the *B neutrino flux were partly due to the strong temperature dependence of this 
flux (see section 4) and to the separate influence of numerous ingredients. This strong 
temperature dependence raised the hope that neutrinos could be a good thermometer 
of the center of the Sun. 

In parallel, helioseismology was being strongly developed on both the theoretical 



and experimental sides ( 


Duvall , 


1979 


Grec, Fossat & Pomerantz, 


1980 


Claverie et al. 


1981 Duvall, 1982 


Christensen-Dalsgaard |1982 


Christensen-Dalsgaard et al. , 1985 



Gough, 1985) with a real potential to properly describe the radiative zone and, maybe. 



the deep core of the Sun. Considering how these two probes (helioseismology and the 
neutrinos) complement each other, S. Turck-Chieze (1988) decided to use the solar sound 
speed obtained from helioseismic instruments in space to go beyond a purely theoretical 



prediction of the neutrino fluxes. Also, she contributed to build the GOLF (Gabriel 



et al. , 1995) instrument, onboard the SoHO spacecraft, dedicated to research on the 



solar core. With her collaborators, she developed more and more accurate standard 
and then seismic solar models. For more details, the reader is refered to the review of 



Turck-Chieze et al. (1993) describing the different aspects of the neutrino puzzle and of 



the solar research. The SoHO spacecraft, launched in December 1995, provided solar 
physicists with a wonderful opportunity to determine the solar sound speed down to 
0.06 (where Rq is the solar radius) with a good accuracy, thus reaching the region of 
neutrino emission (Turck-Chieze et al. |2001b ). In the meantime, most of the physics of 
the standard solar model was improved (see sections 4 and 5). Today, the solar central 
temperature is estimated with a precision level of about 0.1 xlO^ K, much better than 
1 %, thanks to the combined efforts of astrophysicists and the SNO (Sudbury Neutrino 
Observatory) collaboration. 

Since the middle of the nineteen-eighties, the amount of information available 
dramatically increased on both the solar physics and neutrino sides. Different 
neutrino detection techniques were developed to complete the remarkable results of 
the Homestake experiment. Twenty three years after the first result of Davis, a real- 
time experiment began at Kamiokande in Japan, with pure water in an installation 
dedicated to the determination of the proton lifetime. Kamiokande confirmed the solar 
neutrino deficit (Hirata et al. , 1989). Two radiochemical detectors using gallium were 



ready to run at the beginning of the nineteen-nineties, in Gran Sasso (Italy) under 
the responsibility of T. Kirsten (Anselmann et al. 1992; Hampel et al. 1998), and 



in the Balkans under the responsibility of V. Gavrin (Abazov et al. 1991; Gavrin et 



aL) |1992[ [Abdurashitov et al.[ |1994[ |1999a||bl ), thus allowing a comparison of different 
neutrino detection techniques. For the first time, the neutrino flux issued from the 
proton-proton (pp) reaction was measured in addition to the other neutrino sources. 
Strangely enough, the first flux announcements were values just above the predicted pp 
neutrino fluxes for one experiment, and an absence of detection for the other. With an 
increase in the sample size, both experiments showed a global deficit. Then, on April 
1*^ 1996, it was announced that the detector of Superkamiokande started to detect 
solar neutrinos (Fukuda et al. 1998 2000). With Superkamiokande, the statistical 
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accuracy was dramatically improved with more than 5000 neutrinos detected per year, 
coming from the Sun. A dream come true for R. Davis. Superkamiokande was the 
first experiment to unambiguously confirm the existence of neutrino oscillations, with 
the additional measurement of the atmospheric neutrinos. The neutrino oscillation in 
matter as solution to the solar neutrino puzzle was proposed by particle theoreticians 



as early as the first detections by Davis, and has been mainly developed by Wolfenstein 



(1978), and Mikheyev & Smirnov (1986). Therefore, it has been coined the "MSW" 
effect. Finally, the Canadian heavy water detector, SNO, unambiguously proved its 



existence. SNO measures the interaction of different flavors of neutrinos (Ahmad et al. 



2001, 2002 Ahmed et al.l 2004) and allowed, for the first time, a determination of all 



the neutrinos (all flavors) coming from the Sun. The total neutrino flux thus measured 
immediately agreed with the predictions of solar models: standard and seismic models 



of the Sun, whose predictions were close to each other ( 


Bahcall, Pinsonneault & Basu 


2001 


Turck-Chieze et al. 


2001b 


Couvidat, Turck-Chieze & Kosovichev 


2003 


)• 



Undeniably, 2001 was a second very important year for solar neutrino problem. Such 
a wonderful agreement strongly validates the measurements of the pioneer experiments. 
Their authors, R. Davis and M. Koshiba, received a joint Nobel prize in 2002 for their 
seminal contribution to this field. Recently, improvements on solar CNO (Carbon, 
Nitrogen, Oxygen) abundance determination put a renewed pressure on the standard 
solar model neutrino predictions and emphasized the need to recourse to a seismic solar 
model for these predictions. It is clear that solar observations benefit from two probes 
(neutrinos and helioseismology) , and that the Large Mixing Angle (LMA) solution is the 
favored solution to the solar neutrino oscillation ( Bahcall fc Pena-Garay[ 2004 Inoue 



2004 Aharmim et al. , 2005). Helioseismology and particle physics are now in agreement 



(Turck-Chieze et al. , 2004a) at a 10% level, showing the limits of the neutrino-matter 
interaction effects, and the maturity of the different scientific approaches. Borexino 
recently confirmed the impact of the MSW effect by measuring, for the first time, 
the Be and neutrino fluxes with the same detector (Arpesella et al. , 2008). This 
experiment confirms the LMA solution of oscillations, and directly shows the reduction 
of neutrino fluxes at low energy (due to oscillations in vacuum) and at high energy (due 
to oscillations in matter). Physicists are now able to obtain individual measurements 
of the different sources of neutrinos. Borexino and KamLAND (the largest low-energy 
antineutrino detector presently measuring antineutrinos from nuclear reactors) confirm 
the reduction of neutrino flux at low energy due to neutrino oscillations. We will 
probably gain further constraints on CNO abundance in the solar radiative zone in 
a near future, thanks to, among others, the detection of CNO-cycle neutrinos. This 
could be an opportunity to look for magnetic interaction between neutrinos and matter. 
In astrophysics, helioseismology opened the door to the study of extra phenomena not 
yet fully introduced in the standard solar model. For instance, the role of rotation and 
magnetic field in the Sun is beginning to be considered in more detail than what was 
done in the past, which should help bridge the gap in the interpretation of the different 
solar activity data. In the meantime, the seismic solar model remains superior to the 
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standard solar model to explain the different neutrino ffuxes observed (see sections 5, 6, 
and 7). 

The quest for solar neutrinos already produced fantastic results, and is far from 
over. New paths are being explored to answer the fundamental questions of neutrino 
physics: are they Majorana particles, what are their masses? New paths are also being 
explored regarding the Sun and other stars. These two research fields will contribute to 
the development of Neutrino-stellar Astronomy (section 8). 

In the next section we detail the physics of the neutrinos, and its development since 
the "historical time" of the first neutrino flux measurements by R. Davis in 1968. 



3. The weak interaction and the general neutrino properties 

The neutrino is a very peculiar particle: 

- neutrinos are generated in the Sun mainly by the transformation of a proton 
into a neutron, through weak interactions: either /3+ disintegration (with production 
of a positron) or electron capture. The main solar energy source is the interaction 
between two protons which produces a deuteron plus a positron and a neutrino (table 
1). It is remarkable that a neutron (mass of 939.56 MeV/c^) can be generated from 
a proton (mass of 938.28 MeV/c^) despite its smaller mass. The opposite case, the 
(3— disintegration n — )■ p + e + z/g, is conceptually quite natural and free neutrons are 
unstable, with a mean life of 885.7 s (14.7 mn). On the contrary, the long lifetime of the 



proton, estimated from the SuperKamiokande detector to be > 6.6 xlO'^'^ yrs (Nishino 



et al. , 2009), does not immediately favor the reaction: p — )■ + e"*" + i/g- Indeed, it 
is noticeable that a free proton needs an initial velocity (kinetic energy) to produce 
a neutrino, or a reaction between two nuclei yli=A(Z,N) and y42=A(Z-l, N+1) with 
M(Ai) < M{A2)+ rrie. In the Sun, Ai is the sum of two protons and A2 is a deuteron. 
This immediately shows that the cross section of such an interaction will be very small 
and that, consequently, the resulting neutrino energy will be low (< 0.4 MeV). Hence 
the difficulty to detect neutrinos. 

- the absence of neutrino detection by Davis near nuclear reactors, combined to the 
actual detection of solar neutrinos, seemed to favor the idea that Ue (particles emitted 
together with a positron in decay) are different from z/g neutrinos (particles associated 
to electron in /3~ decay, near reactors). This premature conclusion was linked to the fact 
that the neutrino (resp. antineutrino) has a spin 1/2, no electric charge, and a leptonic 
charge +1 (resp. -1). However this simple picture was complicated by the demonstration 
of C. S. Wu (1957) regarding (3 decay: the parity conservation is violated. This followed 



the possible occurrence of such a violation, suggested by Lee & Yang (1956). In the 
minimal standard model of strong and electroweak interactions, SU(3)*SU(2)*U(1), 
the consequences of this violation are the following: the neutrino is a Weyl particule, 
represented by a spinor at two components, and only the neutrinos of left chirality and 
antineutrinos of right chirality exist. The neutrinos can only exist in a specific state of 
polarisation: only left-handed helicity neutrinos are produced (that means that the spin 



Solar neutrinos, helio seismology and the solar internal dynamics 



11 



and the momentum of a neutrino are opposite) and (5~ produce antineutrinos of right 
hehcity (same direction for the spin and momentum). 

- there are 3 famihes of neutrinos associated to the corresponding leptons: electron, 
muon, and tau. The discovery of neutral currents in 1973 showed that all the different 
neutrinos do not experience the same interactions. The different may experience 
neutral interaction mediated by Zq (Z Exchange Process, ZEP), but only v^. can 
experience charged current interaction (W Exchange Process, WEP) mediated by the 



W particle, see Bouchez (2005). 



The idea that neutrinos could oscillate between their different flavors was proposed 



by Gribov & Pontecorvo (1969) soon after the first results of the chlorine experiment. 

, and the chlorine 



Indeed the nuclear reactions in the Sun produce only one flavor, 
detector is sensitive only to this flavor. If there are less neutrinos detected on Earth than 
neutrinos produced in the Sun, this may be due to the transformation of z/e into or 
. One condition for this to happen is that neutrinos are massive. The propagation 
eigenstates are then different from the flavor eigenstates and the phase change during 
the propagation will induce a flavor change. From the neutrino oscillation formalism 
in vacuum and in matter, it was shown how this can be a natural interpretation of the 
solar neutrino experiment results. 



3.1. Vacuum neutrino oscillations 

In the standard model of particle physics, neutrinos are assumed to be purely left handed 
and massless. The standard model is based on arbitrary parameters and is generally 
only considered as a first step toward a more complete theory. In Grand Unified Theories 
(GUTs) neutrinos have a mass and the mass eigenstates z/i , z/2 , and are different 



from the flavor eigenstates v^ 



and Uj- . The total lepton number is conserved, but 
The transformation between the 



not the separate electron, muon, and tau numbers, 
flavor and the mass eigenstates can be written : 

where U is a 3x3 unitary matrix. In the simple case where we consider only two flavors, 
the transformation has only one parameter, 6, which is called the mixing angle, and U 
is written : 







' cos6' 


sin6' 


:: 


-1 


I — sin6' 


cos6' 




When neutrinos propagate in vacuum, the evolution equation takes the following 



form: 



dt 



4E^ 



cos2e 
"sin2^ 



4E^„ 
Am 
4E, 



sin2^ 
-cos2^ 



z/„ 
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where is the neutrino energy, and Am^ = m\ — is the difference of the squared 
masses of vi and z/2 . The probability P(z/c — >z^c ) that a produced at t=0 is still a 
Pe after a propagation time t (or a propagation distance 1) is expressed as: 



P(z/e ^Ve) = I < l^e(t)|l^c(0) > 



1 - -sm^2e 
2 



1 — cos- 



27rl 
17 



(3.1) 



where 1^ = is the vacuum oscillation length. The amplitude of the oscillation 

depends on the mixing angle 6. 

The general solution for three-neutrino flavor is : 



P(z/„ Up) 



5^U^iexp(-iEit) U: 



J]|U^ir|U«i|2 + ReJ]U/3iU^jU;iU„jexp(-27ri-) 



3.2. Neutrino oscillations in matter 



Wolfenstein (1978) observed that the presence of matter modifies the propagation of 



neutrinos because of the effects of coherent forward elastic scattering. As previously 
mentioned, ZEP contributes to the elastic scattering of all neutrinos, whereas WEP 
contributes only to Ue scattering. This implies a difference in the index of refraction for 
Ue and for z/^ or z/^- . The propagation equation is then written : 



^sin2^ 



dt 



^'^-sin2^ 4^cos2^ 




where G is the weak interaction Fermi constant and p is the electron density of the 
medium. 

If the density is constant, the previous system can easily be solved. The probability 
P(i^e -^^e ) has a form very similar to the probability previously obtained in vacuum : 



P(l/e ^ Z/e) = I < Z/e(t)|z/e(0) > 



1 - -sin^2^^ 



1 — COS- 



27rl 



(3.2) 



The mixing angle in matter 9^ and the oscillation length Im depend on the electron 
density of the matter through the relations : 

sin2^ 



cos2e - (W/lo) 
1. 



Vl - 2(W/lo)cos2^ + (W/lo)2 



where lo = 27r/(v^Gp) is a characteristic length of the medium. Mikheyev & Smirnov 



(1986) pointed out that, for a given set of the neutrino oscillation parameters 9^ and 
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Am^ , and for a given neutrino energy , there is a value of the density which induces 
a "resonant" mixing, i.e. sin^26'jn = 1. 

If the density of the matter is not constant, there is no general solution to the 
evolution equation. The mass eigenstates are no longer eigenstates of the Hamiltonian. 



However, if the density has a slowly varying behaviour, Mikheyev & Smirnov (1986) 
showed that the instantaneous mass eigenstates become eigenstates of the Hamiltonian. 
This property, well known in quantum mechanics, is called the adiabatic approximation. 
The crossing of the resonant density may then induce a very important phenomenon: 
the adiabatic transformation of a i^e , flavor eigenstate, into a z/2 , vacuum eigenstate. 
This is called the Mikheyev, Smirnov, and Wolfenstein (MSW) effect. 

This situation may happen in the Sun where the density decreases relatively slowly 
with the radial distance, with an exponential-like behaviour, from about 150g/cm^ at 
the center, where the neutrinos are produced, to almost zero at the surface. In this 
case, a pure z/g can leave the Sun as a pure z/2 , i.e. a pure state of propagation in 
vacuum. Two conditions are necessary for this to happen: a) the density where the 
neutrino is produced must be larger than the resonant density corresponding to its 
energy and to the neutrino oscillation parameters; and b) the adiabatic condition must 
be satisfied. In this case: a) the interval in which there is a flux suppression is large 
and its width decreases with sin^2^ ; and b) the minimum value of the flux is equal to 
sin^^m, which means that the smaller the mixing angle, the larger the Ve flux reduction. 
The reduction factor depends on energy and a consequence is that, depending on the 
oscillation parameters, the observed spectrum may be distorted compared to the 
theoretical one. 

Moreover, when neutrinos reach the detector on Earth at night, there may partly 



be regeneration of the z/g which "disappeared" in the Sun, see Bouchez (1986). Indeed 
if the neutrino arrives on Earth as a pure V2 , this 1^2 is no longer an eigenstate of 
propagation in the Earth and some z/i component may appear depending on the values 
of the parameters. 

3. 3. Characteristics of the different solar neutrino sources 

Neutrinos are produced in the central part of the Sun by several reactions described in 
tables 1 and 2: the pp reaction (with a small contribution, 0.25%, from the weaker pep 
reaction), the electronic capture on ^Be which produces z/bs , and the proton capture 
by ^-Be, which produces excited atoms which decay and generate the so-called z^b • 
The reactions from the CNO cycle involving ^^N and ^^O are of minor importance for 
the luminosity of the Sun but are not negligible for the neutrino flux. The ^^F source, 
which also comes from the CNO cycle, is much less important. A reaction, energetically 
negligible, but which involves ^He and produces neutrinos with an energy up to 19 



MeV (called hep neutrinos) has been identified by Bahcall & Ulrich (1988). The energy 
spectrum (figure 1) is crucial for neutrino detection. This spectrum is determined by 
the kinematics and the energy released by the corresponding nuclear reactions. It was 
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Table 2. The nuclear network: nuclear energy produced by each reaction, maximal 
value of the neutrino energy (E,y max), its mean value deduced from neutrino energy 
spectra (Ej^), for contributions to the total luminosity greater than 0.1%. From Turck- 



Chieze et al.| ( |1988| ). 



Reaction 


nuclear energy 


El, max 


K 


luminosity 




MeV 


MeV 


MeV 




p(p,e+z/)D 


1.442 


0.420 


0.265 


8.26% 


p(pe ,i/)D 


1.442 


1.442 


1.442 





D(p,7) -^He 


5.494 






38.29 % 


S TT TT O \ 4 TT 

•^He[^He,2Y)) He 


12.860 






42.84 % 


^He[a,'y) 'Be 


"1 or' 

1.586 






0.77 % 


7 D„/'„— , ,\7 T A 

\De[e ,u) Li 


0.86/ 0.O/4 




0.86/ 0.3/4 


0. 


^Li(p, a ) ^He 


17.348 






8.07 % 


^Be{p, -f)^B{e+iy) ^Be*{a)'^He 


17.98 


14.02 


6.71 


0. 


i2C(p,7)i3Ar(e+zy)i3c 


4.454 


1.198 


0.707 


0.26% 




7.551 






0.56 % 




10. 05 


1.173 


0.997 


0.60 % 




4.966 






0.33% 


i5Ar(p, 7)160 


12.128 






0. 


i60(p,7)i^F(e+z/)i^O 


3.422 




0.999 


0. 




1.193 






0. 



calculated by Bahcall & Holstein (1986) and updated by Bahcall & Ulrich (1988) for 
other sources. 

All these reactions generate electron neutrinos z/e, but these neutrinos differ by their 
region of emission, their flux, and their energy spectrum (figure 1). The main interest 
of the neutrino lies in its particularly low cross section (about 10~^^ to 10"^^ cm^ in the 
energy range considered: to 14 MeV). This results in a mean free path of about 0.5 
parsec in the solar interior; accordingly, neutrinos might reveal the thermodynamical 
state of the emission region, and detecting different sources of neutrinos might give 
information about different regions of the Sun. 

The regions of emission are determined by the nuclear reaction rates and their 
temperature and density dependence. For the central conditions in the Sun, the three 
major neutrino fluxes can be expressed in the following way ( |Gough|p88| ): 



$pp oc pcXcTc; 



X, 



;i + Xc) 



11.5. 

C ) 



B OC pc- 



(1-X, 



C) rp2A.5 
J- r 



1 + Xc) 



c 



(3.3) 



where Tc is the central temperature, pc is the central density, and Xc is the central mass 
fraction of hydrogen. It should be noticed that these statements are, for all practical 
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purpose, independent of the solar model and only reflect the behaviour of each nuclear 
reaction. 

Figure 2 shows the regions of emission: the ®B neutrinos are produced in the very 
central part of the Sun (the inner 10% in radius or 20 % in mass) because of the 
extreme dependence of this nuclear reaction rate on the temperature; the ^Be neutrinos 
appear in a wider domain, and the pp neutrinos originate from the whole nuclear region. 
The neutrino flux from each source must be integrated over the relevant region of 
emission, which corresponds typically to a region where the temperature decreases by 
20% for neutrinos and more than 30% for pp neutrinos. Thus, the real temperature 
dependence of the neutrino flux, which varies with the solar model, is less steep than 
quoted in Eq. (3.3). Moreover a feedback is always necessary to adjust the luminosity 
of the solar models at present age to the observed value. It turns out that, owing 
to partial compensations due to the variation of the composition and the density in 
different models, the dependence of the solar neutrino flux on the central temperature 
is in fact: 

$ppOC $7^, oc T^; oc T^^ (3.4) 



These expressions were deduced by Bahcall &; Ulrich (1988) from 1000 different 
computations of a solar model with different input parameters covering a large range of 
values. 



4. The classical view of the Sun through the Standard Solar Model 

The time scale required for biological evolution on Earth is considerably greater than 
the Kelvin Helmholtz time of 30 Myr, corresponding to the estimated lifetime of the 
Sun if gravity alone compensates for the loss of luminosity at the surface: this fact was 
at the origin of public disagreements between Darwin and Lord Kelvin at the end of 
the 19th century. At the beginning of the 20th century, it became obvious that the Sun 
should be at least as old as the Earth. This is why Eddington suggested in 1920 that, 
perhaps, a subatomic source of energy may exist in stars. It then became evident that 
thermonuclear fusion had to play a fundamental role. Today, the origin of the Sun's 
stability is clearly established and it is reasonable to treat the long-term evolution of a 
star as a succession of static equilibrium models. 

4-1. The fundamental equations 

The classical picture of the Sun (used to produce "standard solar models" ) was developed 
in the course of the 20th century with the following hypotheses: the star is spherical, 
described by a succession of hydrostatic equilibria, and without effects of rotation and 
magnetic field. We will show in sections 5 and 7 that we are currently expanding 
this picture to better describe the real Sun. The precise knowledge of the general 
characteristics like distance, mass, age, radius, and luminosity, gives the Sun a privileged 
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Figure 2. Spatial emission of the neutrino sources computed with the seismic model. 
In the upper figure are drawn the p-p (plain curve), (dot-dashed curve), ''Be (dashed 
curve), and the pep (dot-dot-dot-dashed curve) neutrinos. In the lower figure, the ^^N 
(plain curve), ^^O (dashed curve), and ^^F (dot-dashed curve) neutrino production are 
shown. For each neutrino type, we have drawn (1/Ft) (dF/dr) where F is the flux in 
s~^, r the fractional radius, and Ft the total flux for this neutrino type (integrated 
over the entire Sun). From Couvidat, Turck-Chieze & Kosovichev (20031. 



position among the stars. These characteristics are summarized in Table 3. Their role 
is crucial in solar modelling because they serve as constraints: for a given model, a 
correct luminosity and radius must be reached at the correct age. In this framework, 
the internal structure of the star is described by four structure equations which are the 
foundations of what is generally called the Standard Solar Model (SSM). 

The first two equations assume hydrostatic equilibrium (each gas shell is balanced 
by the competition between the downward gravitational force and the outward pressure 
gradient) and mass conservation. P, T, and p, are respectively the pressure, temperature 
and density at the position r (radial distance from Sun center), and M(r) is the mass 
enclosed within a sphere of radius r. 

dP M(r)G dM . , , 

— = hHP and — = 47rrV 4.1; 4.2 

dr dr 



Thermal equilibrium is assumed. The energy produced by nuclear reactions (47rr^pe, 
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Table 3. Evolution of the solar fundamental constants and their variability during the 
11 year solar cycle. The heavy element mass fractions have been estimated by |Anders| 
& Grevesse|(|1989|);|Grevesse et all (|1992|); lAsplund et all (|2009|). 



Reference values Allen 



Present values 



Time Variability 



Luminosity 

Age 

Radius 

Seismic radius 
Radius shape 
Heavy element Z 
Present mass loss 



1360.488 (±2.10-^) 



695 990 km 



0.02 then 0.0173 



1367.6 W/m2 - 1361 W/m^ 
4.6 ± 0.02 Gyr 
693710 (min) 
695660 km 
oblateness 6 to 10 km 
0.0134 
2 X lO^i^Mo/yr 



1-4 W/m2 

10-160 km 

6-14 km 
no evidence 



where e is the nuclear energy production rate), balances the energy flux L(r) emerging 
from the sphere of radius r. It includes the energy loss by neutrinos (which is only 
about 0.4 MeV for hydrogen burning, but larger for, say, a supernova explosion). Taking 
into account quasi-static gravitational readjustment and composition variation, a heat 
transfer term TdS is included, where S is the total entropy per gram of the gas and the 
energy conservation yields: 

dL .of TdS\ , , , 

4vrrV e - — (4.3) 



dr \ dt ^ 

The radial temperature gradient depends on the different processes which contribute 
to the energy transport. In a radiative region of a star, the diffusion approximation is 
appropriate, and the relation between temperature gradient and luminosity is: 

dT^-3Kp^rl 

dr AacT^Anr^ ^ ' ' 

When the radiative opacity coefficient k increases too much, like in the Sun, or when the 
energy production is very high (in the internal part of stars with mass > I.5M0), the 
radiative gradient increases so much that matter becomes convectively unstable. The 
resulting temperature gradient is then nearly adiabatic: 

dr \dr ) T2 P dr ^ ^ 

where and the other adiabatic exponents, are defined by: 

PV^^ = const ; p^-^^T^^ = const ; TV^^'^ = const (4.5) 

In a monoatomic ideal gas, Fi = r2 = Fa = In stars, these three quantities are not 
equal, especially because of the partial ionization of different elements. Moreover, these 
three quantities drop from 5/3 to 4/3 when a pure ideal gas is being replaced by pure 
radiation. The value F3 < 4/3 is, for example, at the origin of pulsation of stars. Fi 
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determines the dynamic instability, while r2 governs the convective instability, and Fa 
governs the regime of pulsation instability. 

Boundary conditions (typically for pressure and temperature at the stellar surface 
and for mass and luminosity at the stellar center) are included to solve these equations 
and to follow the solar structure evolution. 



4.1.1. The solar equation of state The total pressure P in a star is the sum of radiation 
pressure and gas pressure: 

Pn = laT^;Pa=^kT=P^ (4.6) 
3 V fi{r) 

where a is the radiation density constant and Vis the volume of the gas. The radiation 
pressure Pr is negligible in the solar center (1/1000 of the gas pressure Pq) but not near 
the photosphere (.6Pg)- When temperature is high and density low enough, interactions 
between the particles are negligible and the gas can be reasonably approximated by an 
ideal gas (equation 4.6). This approximation is quite good in the central part of the 
Sun for ions. However, the large range of solar temperatures (5 x 10^ to 10^ K) and 
densities (10~^^ to 10^ g/cm^) requires a more detailed description of the equation of 



state (Rogers, Swenson & Iglesias, 1996 Rogers & Nayfonov, 2002). Since the star is 



constituted of a mixture of different chemical species, gas pressure and mean molecular 
weight fi are written as: 

P = y" —pkT 1^-^ = V ni(mn/mi) (4.7) 

i i 

In these expressions, is the number of free particles: Z electrons + 1 nucleus for 
each atom of atomic number Z. More generally, this number depends on the degree of 
ionization of the species i considered. In astrophysics, X, Y, Z represent the mass fraction 
numbers (normalized to 1) of, respectively, hydrogen, helium and other elements (heavier 
than helium). For a completely ionized gas, /i is expressed as = 2X + (3/4)y + 0.5Z 
where it is assumed that the proton number equals the neutron number for heavy 
elements. The solar value /x increases outward due to the partial ionization of elements 
and the existence of molecules. As shown by their high ionization potential. He and 
Ne are difficult to ionize. Central solar density is sufficiently high to cause partial 
degeneracy of the electrons. Therefore, electron number rie and electron pressure Pg 
must be expressed as functions of the degeneracy parameter r): 

Pe = ^-I^{2mkTf'F,„{-r^) -n, = ^(2mkTf'^Fy^{-,^) (4.8) 

where P1/2 and P3/2 are the Fermi-Dirac functions. In the case of the Sun, the effect of 
degeneracy is small (r/ of the order of -1) but, if the star continues to evolve, it reaches 
a point where this effect becomes dominant. This leads to an increase in the pressure 
and a quasi independence of the pressure on the temperature. 
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4.1.2. The thermonuclear source of energy The calculation of the nuclear energy rate 
e requires a good knowledge of the reaction rates that take place in the solar interior. 



Gamow (1928) was the first to show that hydrogen is the most abundant element in 



stars, and is the first element to be converted because it has the lowest Coulomb barrier 



(energy barrier that two nuclei need to overcome). Von Weizsacker (1938) and Bethe 



(1938) clearly showed that two different sets of reactions, the pp chain and the CNO 
cycles, could provide the amount of energy explaining the present solar luminosity. In 
the solar case, the pp chain is the most energetic one (e oc T^), but when temperature 
increases (at the end of hydrogen burning), the CNO cycles become the most efficient 
(e oc or T^^). This is why low mass stars (< I.SMq), driven by pp reactions, have a 
slower evolution than more massive stars, driven by CNO cycles. 

A typical nuclear reaction between two species a and X can be described by: 

a + X -^b + Y 

The associated gain or loss of energy Q is given by the energy balance Q = AE = 
(Mx + Ma — My — Mb)c^. If cr(t>) denotes the probability that a projectile a collides 
with a fixed nucleus X with relative speed v, the total number of reactions per cm^ and 
per second is given by: r = NaNxV(T{v) where Na and Nx are the number densities 
of the species a and X. Since matter in the stellar interior is in local thermodynamical 
equilibrium, the relative velocities are given by a Maxwellian distribution f{v), and the 
total reaction rate is expressed as: 

r = NaNx / va{v)f{v)dv = NaNx < av > (4.9) 



The total energy generated per unit mass and per unit time is: e = ^ where p is the 
gas density. At typical stellar temperatures, the average thermal energy of a particle 
(considering a Maxwell-Boltzmann distribution of velocities) is several keV, which is 
about 10~^ times smaller than the Coulomb barrier (~ 1 MeV). Therefore, from a 
classical point of view, no reaction could take place, and the pp interaction could only 
occur when the internal temperature reached 10^*^ K (corresponding to the Coulomb 
barrier). However, at that point, the ensuing nuclear reactions would be so strong 
that the star would experience a catastrophic explosion. Of course, nuclear reactions 
could also occur in the high energy tail of the Maxwell-Boltzmann distribution, but 
the number of participating particles would be very small. This problem was solved by 



Gamow (1928) when he showed that, in quantum mechanics, there is a small but finite 
probability (1.7xl0~^ at 5 keV) that the proton penetrates the Coulomb barrier. This 
tunneling probability is proportional to: 

exp[-27raZiZ2(|^)^/^] = exp (-^) (4.10) 

where is the charge of the species i, and a = e^/hc. The generation of energy is 
expressed as: 

iv.ivxo, 8 ,1,2,, „ -3/2 r -E b 



^J"m-"'l Smexp(--^M£ (4.11) 
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where S(E) gives the smooth variation of the cross section with energy E. In fact, the 
cross section for charged particle-induced nuclear reaction drops rapidly for energies 
below the Coulomb barrier and extrapolation is often needed to determine the factor 
(j{E) = exp {—27iri) S (E) . The competition between these two trends, decrease of 
the Maxwellian distribution with E and the increase of the penetration effect, leads 
to an increase in the cross section in a rather small range of energies in the vicinity of 
Eo = (6A;T/2)^/^, known as the Gamow peak. E^ increases with the charge of the particle 
and is typically between 10 to 50 keV for hydrogen burning at a central temperature of 
15 million K (1.3 keV). Nevertheless, to calculate the production of energy of (4.14), all 
coefficients S(E) must be determined from laboratory measurements for the reactions 
involved. In the case of non-resonant reaction rates in a non-degenerate, non-relativistic 
gas, the reaction rate r may be expressed as: 



S 



12) 



■ 2 xi/2 ^-^0 e 



-3^0 
kT 



'To' 



(4.12) 



where AEq is the full width at 1/e of the maximum value of the Gamow peak. The 
effective value 5*6// of the cross section factor, S(E), is given by: 



5(0) 



5kT S'(0) 



^ 35, \ 1^(0) 



MeVbarn (4.13). 



36^0 ' .5(0) V " ' 36" J ' 2 5(0) V ' 36" 
The extrapolations of «S'(0), >S"(0), S"{0) from laboratory measurements were performed 



by Fowler, Caughlan G R & Zimmerman (1967), (1975), (1983) for many thermonuclear 
reactions. They expressed the product Na < av > as: 



Na < av >= CiT-'/'exp[C2T-^/=^ - (Tg/To 



(4.14) 



'1 



C^tI'^ + + C5T9 



CqT^'^ + C^Tl''^) cm=^sec-^(mole)-^ 



where all the coefficients Cj depend on Za, Zx-, and 5'(0) or 5"(0) or 5"'(0) and 
Tg = T/10^. Most of the cross sections relevant to the solar interior were measured and 
extrapolated in the twentieth century. Their recommended values have been commonly 
defined in [Adelberger et "aL (1998). A few specific cross-sections have recently been 
improved (see section 4.3). 



J^..2. Energy transport 

Two processes dominate the transport of energy in a solar-type star: radiation and 
convection. In the Sun, conductivity plays no role whatsoever (it has an effect in higher- 
density stars and in supernova cores) and the energy transported by the neutrinos is 
very small. In the central part, radiation is the most efficient transport process. 



4-2.1. Radiative transport In the solar interior, the average energy per photon is in 
the keV range, which is characteristic of X rays, whereas the light escaping from the 
surface is dominated by the visible spectrum, corresponding to an energy 10^ times 
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lower. The cause of this sharp decrease in photon energy is the couphng between 
radiation and matter. The effect of the interaction of photons with electrons, atoms, ions 
and molecules must be evaluated and is summarized by what is called the "Rosseland 
mean opacity" k, which is a function of temperature, density, and composition. The 
equation of radiative transfer adapted to stellar interiors is given by equation (4.4a). 
The radiation field is nearly isotropic and Planckian. The opacity of the gas really 
drives the rate at which energy escapes, and plays a crucial role in determining the 
age of a star and its central temperature. The sources of opacity are the following: 
(1) electron scattering (Thompson and Compton scattering): scattering of the photon 
on free electrons; (2) bound-bound transitions: the photon triggers a change of energy 
level of a bound electron; (3) bound-free (photoionization) transitions: a bound electron 
becomes free after interacting with the photon (or the opposite); (4) free- free transitions 
(inverse bremsstrahlung) : a free electron stays free after the interaction with the photon; 
and (5) molecular transitions at the photospheric level. 

The computation of all these processes necessitates detailed calculations of complex 
atomic physics. Frequency-dependent opacities are needed for detailed stellar models: 
i^toti^) = where i represents the different processes. It is noticeable that 

the Rosseland mean opacity k is not obtained by averaging the frequency-dependent 
opacities, but by averaging their reciprocals, which are the radiative conductivities (the 
weight function /(z/) is the temperature derivative of the Planck function): 



(4.15) 



The major difficulty with such calculations arises from the fact that a single source 
of opacity may exceed all the others. For example the bound-bound processes may 
introduce cross sections which are 1 or 2 orders of magnitude greater than the competing 
contributions. 

Therefore, even a trace element may produce a non negligible contribution to the 
Rosseland mean opacity. In the case of the Sun or solar like stars, the detailed knowledge 
of radiative processes is crucial. Indeed, the contribution of heavy elements (everything 
heavier than helium, representing less than 2% in mass fraction and 0.14% in number 
fraction) to the equation of state does not exceed 1 or 2%, but their contribution 
to the opacity is about 30% in the solar central region, and more than 70% in the 
intermediate region. This is why the iron contribution reaches 20 to 30 % of the total 



opacity (see Turck-Chieze et al. (1993, 2009b) for more details). In the latter reference 



we discuss how these estimates can be evaluated by measuring the opacities, in some 
specific cases, with high-energy lasers. All the radiation-matter interaction processes 
must be calculated for all the chemical elements. Moreover, when peaks appear near 
the maximum of the Rosseland weighting scheme, the bound processes play an even 
more important role. At low energies, more levels contribute to the opacities, and the 
number of peaks greatly increases. Clearly, since bound-bound and free-bound processes 
vary as Z^, instead of for free processes, the heavy elements play an important role 
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Figure 3. Total solar opacity coefRcient as a function of the radial distance from Sun 
center. This figure shows the impact of the different elements, mainly through bound 
processes, on the increase of the opacity. The most important elements are clearly 
visible : C, N, O (at the limit of convective instability) and at low temperature He 
and H. The appearance of neutral atoms or even molecules produces the decrease of 
the opacity near the surface. From |Turck-Chieze et al. (19931. 



far exceeding what could be expected from their low abundances. Moreover, due to 
the non-additive character of the opacity coefficients, the spectrum must be calculated 
for the correct chemical element mixture before deriving the Rosseland mean value. 
This leads to extensive calculations which must be done carefully for different mixtures. 
Figure 3 illustrates the specific contribution of different elements when they become 
partially ionized and when bound processes start to significantly impact the opacities: 
first (in term of distance from Sun center) Fe, then Ne, Si, O, N, C... When oxygen 
(the most abundant heavy element) becomes partially recombined, the temperature 
gradient greatly increases. Further away from Sun center, nitrogen and carbon exhibit 
a similar behaviour and the resulting increase in k becomes so significant that the 
radiative temperature gradient exceeds the adiabatic one: 



dT 




dT 






> 






dr 


rad 


dr 


ad 



(4.16) 



This is the criterion for the onset of convective instability, refered to as the Schwarzschild 
criterion. Therefore in the Sun, there is a large convective solar region below the 
photosphere. In this region conditions are satisfied for first (in term of radial distance 
from solar center) oxygen, then nitrogen, carbon, and finally helium and hydrogen, to 
significantly impact the opacities through bound-bound contributions. However, closer 
to the solar surface, at about O.995i?0, hydrogen and helium become neutral, and the 



Solar neutrinos, helio seismology and the solar internal dynamics 



23 



opacity decreases, prompting the radiative flux to become significant again: we enter in 
a superadiabatic regime. 

4.2.2. Convective transport In any case, the energy produced at Sun center must 
escape from the star. When the opacity increases too much, a very steep temperature 
gradient is required to maintain the energy flow, leading to an unstable situation. When 
the instabihty criterion is fulfilled, the real temperature gradient exceeds the adiabatic 
gradient and the difference is called the superadiabatic gradient: 

MT dT AT^ , AT^ - 1 T rfP 

^ " ^ - ^^^^^ ^^^-^ = I-Pl^ ^^-'^^ 



Convection is a very efficient means of energy transport. In fact, near the base of the 
convective envelope, the radiative flux is so small that we can ignore the superadiabatic 
gradient and use the adiabatic one as the temperature gradient (this is also the case 
in the central part of massive stars). Unfortunately, this is not the case in the outer 
layers of the Sun. In ID solar models, the convection is usually treated by the "mixing 



length approximation" (Vitense, 1953 Bohm, 1958) which replaces the real situation 



of "plumes", convective eddies of different sizes, etc... with an average situation where 
each convective element travels a distance A, before mixing with the surrounding matter. 
This distance is generally scaled to the pressure scale height [|] by a parameter a called 
the mixing length parameter: 

,dlnP~^ 

A = aX = a(^—) 4.18 
dr 

The relative amount of energy carried by convection and by radiation must be estimated. 
The temperature of a convective element and the average temperature of the surrounding 
matter are very different in the central and superficial solar regions. In the central region, 
the temperature gradient is about 10~^K/cm with a mean temperature of about 10^ K 
and a convective time of about 1 to 100 days. Therefore if a convective instability exists 
in the core (for M> 1.5 Mq), it is an extremely efficient process of energy transport 
and the temperature gradient can be approximated by the adiabatic one. In the outer 
layers, the temperature variation between shells is of the order of 100 to 1000 K (to 
be compared with a gas temperature of 5000 to 10000 K), the convective time is much 
shorter (1 to 100 minutes), and we really need a specific treatment of the convection. In 
that case, due to partial ionization of hydrogen or helium, the second adiabatic exponent 
r2 drops to a value close to one and the adiabatic gradient becomes relatively small. 
This favors convective transfer but the radiative flux remains important and cannot be 
neglected. Convective models have been greatly improved by the extension of the mixing 
length approach: some models take into account different Kolmogorov energy cascades 



(Canuto, Goldman & Mazitelli, 1996). Moreover, 3D simulations of such convective 



effects are now at a very advanced stage and reproduce realistically the surface and 



atmospheric turbulence (Nordlund, Stein & Asplund, 2009). 



The pressure scale height is the distance corresponding to a decrease in pressure by a factor e ~ 2.718. 
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The value of the mixing length parameter is still calibrated in solar models, so 
that these models match the solar radius at the solar age. Therefore, in this specific 
case the small superadiabatic region just below the photosphere can be considered 
as correctly treated. The use of this solar value for other stars with rather different 
surface temperatures was considered more uncertain due to the fact that the mixing 
length parameter is largely dependent on the opacity in the superadiabatic part of the 
convection zone. It is now possible to test this approach by comparing precise radius 
measurements of red giants to models using different theories of convection (Piau et 
al. 2010). Significant progress will also be achieved with stellar seismology (especially 
with the COROT and KEPLER missions), which allows the derivation of the convective 



depth for various solar-like stars (Ballot, Turck-Chieze & Piau, 2004). 



4-3. Improvements in the physics of the SSM. Evolution of the neutrino flux 

Radial temperature, density, pressure and composition of the present Sun (and 
consequently the emitted neutrino fluxes) are obtained by solving the equations (4.1) to 
(4.4), at each time step, as a succession of hydrostatic equilibria. The SSM is the most 
economical way of describing the Sun: it assumes that there is no important effect of 
rotation and magnetic field. In this framework, only 3 observed variables are needed: the 
luminosity of the Sun at an age of 4.55 Gyr (from the onset of hydrogen burning, or 4.6 
Gyr since the solar formation), the solar radius, and the detailed element composition for 
nuclei greater than helium (determined at photospheric level and compared to meteoritic 
compositions). The only free variable is the initial helium content which is adjusted 
in order to reproduce these observations at the present time. Solving the equations 
(4.1) to (4.4) requires the knowledge of the physical inputs previously described, to 
derive P(p, T, Xi)(r), e(/9, T, Xi)(r) and k(p, T, Xi)(r) at each time step. Knowledge of 
these physical ingredients has been significantly improved during the last two decades. 

Furthermore, several of these improvements have been accomplished since 1988, 
the year when the European Space Agency decided to favor helioseismology from space. 
Table 4 summarizes the physical processes introduced in the solar models which can be 
estimated thanks to one or two observable quantities, and the location in the Sun where 
a specific element reaches its peak impact. 

4-3.1. The detailed solar composition: Xi(t,r) The detailed composition of the plasma 
is a key ingredient of the SSM. It enters in the calculation of radial pressure, density, 
and opacity coefficients, as previously shown. It determines the duration of the proton 
burning phase which is largely influenced by the transparency of the star. The knowledge 
of the composition is required at various stages of the theoretical approach. The 
elements from hydrogen to oxygen determine the evolution of nuclear burning; the 
initial composition also influences the mean molecular weight, which plays a crucial role 
in the determination of the pressure. Consequently, the neutrino flux is sensitive to 



the metal composition, as shown by Bahcall & Shaviv (1968). The solar abundances 
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Table 4. Location in the Sun where a specific element can be used to shed light on 
a physical process thanks to a specific observable quantity. 

reference 



radius 

R© 


1 1 

element 


physical process 


1 11 

observables 


0.98 


^He 


microscopic diffusion 




0.71 


160 


transition radiation/convection 


C^, K 


0.70 


^Li 


nuclear process, turbulence 


(? , rotation 






role of magnetic field 


c^, rotation 


0.57 


^Be 


nuclear process, turbulence 




0.25 


3He 


presence of mixing ? (no) 




0.05 


^Be 


presence of mixing ? (no) 


c^, neutrinos 


O.- 




central temperature 


c^, neutrinos 


O.I 









Vorontsov, Baturin &: Pamiatnykh ( 1991 ) 



Christensen-Dalsgaard &: Berthomieu 



Brun, Turck-Chieze & Zahn (1999) 



not solved 



(1991) 



Brun, Turck-Chieze fc Zahn (1999) 



Turck-Chieze et al. 



Turck-Chieze et al. 



(2001a) 



(2001a) 



iTurck-Chieze et al.l (|20oiap ;(2004) 



were poorly known at the time Bahcall & Shaviv (1968) was written, explaining why 



the SSM was calculated with different hypotheses leading to a difference by a factor of 
two in the neutrino flux predictions. 

To determine the initial composition of the Sun, three sources are used: the Earth, 
the solar photosphere, and the meteorites. Each source raises some specific issues. The 
Earth has lost a large fraction of its volatile elements and there is a chemical fractionation 
of the various elements in the different terrestrial layers. On the other hand, the Earth 
is very useful in most cases to determine the isotopic ratios which are not influenced by 
the chemical fractionation and can be very accurately determined. Two other sources 
of information are currently favored: the direct spectroscopic observation of the solar 
photosphere and the chemical analysis of CI carbonaceous chondrite meteorites. Both 
are assumed to characterize the initial composition of the protosolar nebula. Therefore 
their agreement, or lack thereof, is extremely interesting. They inform us on the thermal 
and chemical history of the early phases of the formation of the solar system: precise 
determination of the possibly biased meteoritic abundances has to be compared with the 
reliable but often delicate photospheric abundances that depend on the knowledge of 
atomic physics of the solar atmosphere. The convective zone does not penetrate deeply 
enough to alter the surface composition, which is why the photospheric abundances have 
been considered for a long time to be the initial abundances of the Sun for most of the 
elements. In reality, a slow microscopic diffusion leads to a decrease with time in heavy 
element abundances, relative to hydrogen at the surface, of typically 10% in 4.6 Gyr 



(Michaud & Proffitt, 1993). This constraint on the solar models is added to obtain the 
present photospheric abundances at the age of the Sun. Moreover, some elements are 
burned at very low temperature, such as deuterium, lithium, and beryllium. For these 
elements, the photospheric values are affected by the Sun's evolution. 

After years of studies and precise measurements of atomic oscillator strengths 



(Grevesse et al. , 1992), the photospheric determination of element abundances is very 



much in agreement with the meteoritic one, especially for refractory elements (Table 
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Table 5. Abundances of the elements used in a solar model derived from the solar 



Photosphere (from Anders fc Grevesse| ( |1989[ ): AG, from |Asplund et al. ( 2005 ): AGS 
and from Asplund et al. ( 20091) and compared to meteoritic values from Lodders ( 2003 1 . 



They are given in fraction number normalised such that log N// = 12. The main 
changes are a strong reduction for the C, N, and O abundances (see text). Indirect 
solar estimates are mentioned in parentheses. 



EleniGnts 


Pliotospliere AG 


IVIeteoritGS 


Photosphere AGS 


Asplund et al. (2009) 




12.00 




12.00 


12.00 




[10 99 + 0351 


1.2900 


no 93 ± Oil 


\W 93 ± Oil 




1.16 ± 0.1 


3 25 ± 06 


1.05 ± 0.1 


1.05 ± 0.1 


^Be 


1.15 ± 0.1 


1.38 ± 0.08 


1.38 ± 0.09 


1.38 ± 0.09 




8 56 ± 04 


7 40 ± 06 


8 39 ± 05 


8 43 ± 05 




8 05 ± 04 


6 25 ± 07 


7 78 ± 06 

t • I \J —i— \J •\J\J 


7 83 ± 05 


^0 


8 93 ± 035 


1 38 ± 08 


8 66 ± 09 


8 69 ± 05 




[s 09 ± 11 




f7 84 ± 061 


17 93 ±0 11 


iijVa 


6 33 ± 03 


6 27 ± 03 


6 17 ± 04 


6 24 ± 04 




7 58 ± 05 


7 53 ± 03 


7 53 ± 09 


7 60 ± 09 




6.47 ±0.07 


6.43 ±0.02 


6.37 ±0.06 


6.45 ± 0.06 




7.55 ±0.05 


7.51 ±0.02 


7.51 ±0.04 


7.51 ±0.03 


15p 


5.45 ±0.04 


5.40 ± 0.04 


5.36 ±0.04 


5.41 ±0.03 


16^ 


7.21 ±0.06 


7.16 ±0.04 


7.14 ±0.05 


7.12 ±0.05 


^^Cl 


5.5 ±0.3 


5.23 ±0.06 


5.50 ±0.3 


5.50 ±0.3 




[6.56 ±0.1] 




[6.18 ±0.08] 


[6.40 ± 0.08] 


19^ 


5.12 ±0.13 


5.06 ±0.05 


5.08 ±0.07 


5.03 ±0.07 




6.36 ±0.02 


6.29 ±0.03 


6.31 ±0.04 


6.34 ±0.04 




4.99 ±0.02 


4.89 ±0.03 


4.90 ± 0.06 


4.95 ±0.06 




5.67 ±0.03 


5.63 ±0.05 


5.64 ±0.10 


5.64 ±0.10 




5.39 ±0.03 


5.47 ±0.03 


5.39 ±0.03 


5.43 ±0.03 




7.67 ±0.03 


7.45 ± 0.03 


7.45 ± 0.05 


7.50 ± 0.04 




6.25 ±0.04 


6.19 ±0.03 


6.23 ±0.04 


6.22 ±0.04 



5). This agreement is often better than 5 % even if noticeable differences still persist 
on elements like CI, Mn, Fe, Ga, and Ge. Following that compilation, iron received 
specific attention because its photospheric abundance Fehigh = Fe/H = 4.68 ± 0.33 10~^ 
was 30% greater than the meteoretic one. However, Holweger, Heise & Kock (1990) 
based on ionized Fe (95%) derived an abundance much closer to the meteoritic value: 
Feiow = Fe/H = 3.24 ± 0.075 10^^. The consequences of this improvement on solar 
neutrino fluxes have been discussed in Turck-Chieze & Lopes ( 1993[ ) (see their table 4). 

The lack of information on H, C, N, and O in meteorites (because these elements 
do not enter in the composition of condensable solids) raises the problem of the 
normalization of the meteoritic abundances relative to hydrogen. This normalisation 
is in fact performed using meteoritic and photospheric determination of the refractory 
elements (Mg, Al, Si). For volatile elements the deviation is smaller than 20 %. 
Therefore it is hoped that, in the case of very heavy elements (A> 56) for which 
photospheric data are not reliable, CI meteoritic abundances are representative of the 
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solar properties. These elements are not used in the calculation of the SSM or of seismic 
solar models. Therefore, we do not list them in table 5 but their fraction number appears 
in Asplund et al. (2009). It is estimated that their influence is very small, but their 
impact on the opacity coefficients must be verified. 

The C, N, and O elements are the main contributors to the "heavy elements" 
category. These major volatile elements can only be studied in the solar photosphere. 
Their abundances relative to hydrogen are extracted from molecules like CO, CH, OH 
and NH, or from neutral lines. Their abundance determination requires the knowledge 
of accurate oscillator strengths. Several revisions were published in the last decade 



(Grevesse et al. , 1990 1991 Biemont et al. 1991), and a recent decrease by 30% in 



the abundances of these elements has been advocated by several authors (Holweger 



2001 Asplund et al. , 2004): such a decrease results, among others, from the use of 
better spectral lines. As a consequence of these successive updates in the heavy element 
abundances, the solar Z value decreased from Zq = 0.02 to the present value of Zq= 



0.013, in 20 years. Many papers focused on this update (Asplund et al. , 2005 Caffau et 



al. 


2008 


al. 


2009 



2009 ), and the decrease is summarized in the last column of table 5 (Asplund et 



The consequences of these updates on the B neutrino flux prediction is shown 



in table 6, as well as on the sound speed in the radiative zone (Turck-Chieze et al. 2004a 



Guzik, Watson & Cox 2005), rapidly followed by Bahcall et al. (2005). This update 



also started a discussion on the ability to check solar composition by helioseismology 



(Basu & Antia 2004). Today this promising idea remains difficult to implement, due 
to the variation between the existing equation of state (for the external layers) and an 
insufficient radial accuracy on the sound speed in the radiative zone. The detection 
of gravity modes may lead to some progress and to a better constraint on the opacity 
coefficients. Still, the most promising way to determine the inner composition remains 
the use of both helioseismology and the CNO neutrino detection. New MHD calculations 
put in evidence the role of the magnetic field on the broadening of the lines and could 



lead to a slight increase of the C,N, O, Fe photospheric lines (Fabbian, Khomenko 



Moreno-Insertis & Nordlhund, 2010) 



4-3.2. Solar composition and our Galaxy evolution An important result of solar 
modelling is the determination of the pre-solar helium abundance, derived by forcing 
the solar model to reach the present luminosity at the present age. This value is also 
reported in table 4 and is compared to nearby HII interstellar medium (interstellar 
gas of temperature around lO^i^ where hydrogen is ionised) and to the composition of 
hot stars to determine whether the Sun is typical of our neighbourhood or peculiarly 



rich in heavy elements (Peimbert, Storey & Torres- Peimbert 1993). Helium abundance 



cannot be directly measured in the photosphere due to its low temperature, but seismic 
derivation of the Fi profile in the region where helium is partially ionized allowed a 
determination of its abundance immediately below the photosphere. The photospheric 
value of Y = 0.249 obtained by Vorontsov, Baturin & Pamiatnykh (1991) puts a strong 
constraint on the gravitational diffusion of elements during the life of the star (see 
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Table 6. Evolution with time of the SSM or seismic predictions of the neutrino 
flux in fO^cm^^s^^. Added are the central temperature Tc in fO^X, the initial helium 
abundance Y in mass fraction, and a specific problem that was solved. From |Turck-| 



Chieze et al. (2010a 



flux 


Tc 


Y initial 


problem solved 


reference 


3.8 ± 1.1 


15.6 


0.276 


CNO opacity, ^Be(p,7) 


TC88 


4.4 ± 1.1 


15.43 


0. 271 


-30% Fe abundance, screening 


TCL93 


4.82 


15.67 


0.273 


microscopic diffusion 


BTCM98 


4.82 


15.71 


0.272 


turbulence tachocline 


BTCZ99 


4.98 ± 0.73 


15.74 


0.276 


seimic model 


TCOlb 


5.07 ± 0.76 


15.75 


0.277 


seismic model, magnetic field 


Cou2003 


3.98 ± 1.1 


15.54 


0.262 


-30 % CNO composition 


TC2004 


5.31 ± 0.6 


15.75 


0.277 


seismic model+ '^Be and ^^N(p,7) 


TC2004 


4.21 ± 1.2 


15.51 


0.262 


SSM Asplund 2009 


TC2010 



below). This phenomenon of diffusion explains the difference between the helioseismic 
value and the initial value derived from a solar model (see its evolution in table 6). 
The solar model compatible with seismic data favors an initial Y = 0.277. This value, 
combined with the recent 0/X determination, suggests that the Sun was not born in an 



enriched environment resulting from a supernova explosion (Holweger, 2001; Asplund 



et al. , 2004), contrary to what was previously thought. It is also the conclusion of 



Peimbert et al. (2007). A new galactic helium abundance law can be derived: He/H= 
0.075 -|- 44.6 0/H in fraction number. This law is now compatible with all the indicators 



(Turck-Chieze et al. 2004a) called TC2004. This conclusion is also supported by other 



studies on radionucleides (Gounelle et al. , 2001 2003; Gounelle & Melbom et al. , 2008) 



Ji^.S.S. The interaction between photons and matter K,{Xi, p,T) : Photons escape from 
the nuclear region and interact with matter through the phenomena described in section 
4.2. The Rosseland opacity coefficient, /t(Xj,p, T), has a strong influence on the 
central temperature. The opacity calculations require an accurate determination of 
the composition, especially the helium, carbon, oxygen, and iron composition, and 
a detailed calculation of the photon-matter interaction, obtained from a thorough 
knowledge in atomic physics. Astrophysicists commonly use tabulated values obtained 
for specific compositions. Most of the solar models use, for temperatures larger 



than 6000 K, the opacity compilation of Iglesias & Rogers (1996) who considered 21 
elements: H, He, C, N, O, Ne, Na, Mg, Al, Si, P, S, CI, Ar, K, Ca, Ti, Cr, Mn, Fe, 



and Ni. At lower temperatures, the compilation by Alexander & Ferguson (1994) is 



preferred. Recent opacity tables were recalculated based on the new abundances (see 



http: / /webs. Wichita, edu/physics/opacity/ ). 

It is noticeable that the successive composition updates (Fe, C, N, and O) have 
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systematically deteriorated the agreement between SSM prediction of the sound speed 



in the radiative zone and the "observed" (helioseismic) sound speed (Turck-Chieze & 



Lopes, 1993 Turck-Chieze et al. 2004a). It is the case also for the location of the base 
of the convective zone and for the photospheric helium. Three directions of investigation 
appear today: the magnetic effect on the photospheric lines, the detailed understanding 
of the photon interaction, some effort is under way to validate these complex calculations 



Bailey et al. , 


2007 


2009 


Loisel et al. 2009 


Turck- 



Chieze et al. 2009b 2011a), the third is to go beyond the SSM (section 7). 



4.3.4- Temporal evolution of the composition: Xi{p,T) Heavy elements diffuse toward 
the center during the solar life. This microscopic diffusion, mainly gravitational settling 



(GS) in the solar case, is a very slow process (Proffitt &; Michaud, 1991). However, it is 
a crucial phenomenon that must be introduced in the equation describing the temporal 
evolution of the composition: 



dt dt niiri dm, 



■ VjXi 



(4.19) 



dt dt nucl dm 
Microscopic diffusion explains the observed stellar photospheric abundances, and its 



in solar models ( 


Christensen-Dalsgaard, Proffitt & Thompson 


1993 


Thoul 


Berthomieu, Provost & Morel 


1997; 


Brun, Turck-Chieze & Morel 1998) 



improved the agreement between theoretical and observed sound speed (Fig. 6a). In fact 
we introduced two terms in equation (4.19) in addition to the nuclear term: the first one 
describes the migration of elements relatively to hydrogen (microscopic diffusion) and 
the second is a turbulence term in the region of transition from the radiative to convective 



energy transport, which partly inhibits this microscopic diffusion (Brun, Turck-Chieze 



fc Zahn||1999[ ). 

The microscopic diffusion of the elements modifies the composition along the radial 
profile and reduces the hydrogen burning lifetime by almost 1 Gyr. The photospheric 
composition (relative to hydrogen) is reduced by about 12 % during this lifetime. The 
introduction of such a process increased the neutrino flux estimates for the chlorine and 
water detectors by about 20% (table 6). The introduction of such a dynamical effect 
helps to better understand the photospheric abundances of lithium and beryllium, which 
burn at respectively 2.5 x lO^K and 3.5 x lO^K (see also section 8). Surprisingly, this 
change almost compensates for the impact of the change of composition on sound speed, 
temperature, and neutrino fluxes. 



4.3.5. A detailed description of the nuclear interaction in the solar plasma The 
knowledge of the reaction rates listed in Table 1 (except for the very weak pp reaction 
rate) are based on measured cross sections in laboratory. Exceptionally, {^He,^ He) 



1998) 



has been measured down to the astrophysical range of energy (Junker et al. 
Updates were obtained on the (jBe.p) cross section (Junghans et al. , 2003) and on 
(^^iV,p), which was reduced by a factor 2 (Formicola et al. , 2004) compared to previous 
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estimates. However, the reaction rates described by equation (4.12) are modified by the 
solar plasma. The velocities of the reactants are indeed perturbed by the presence of free 
electrons and ion clouds in their vicinity. A correction, called "the screening factor" , is 



needed and was further discussed in the framework of the neutrino problem (Dzitko et 



al. 


1995 


Gruzinov 


1998 



the pure weak plasma described by the Debye theory. The current solar observations 
are compatible with the notion of intermediate plasma and any larger screening effect 
seems unwarranted. The pp reaction rate is so small that it is only known theoretically. 
As the sound speed is extremely sensitive to this reaction rate, it is now possible to 



improve its determination through seismology (Turck-Chieze et al. , 2001b) and also 
through the neutrino flux measurements. Experimental efforts are pursued to improve 
the nuclear reaction rate predictions, and a new compilation, including a discussion 



on the various sources of uncertainty, was published this year (Adelberger et al. , 2011). 
Laser measurements on LM J or NIF represent a new challenge for the future in producing 
plasma in stellar conditions, this new technique has the advantage to measure directly 
the nuclear reaction rates. 

The possibility that physical mechanisms such as diffusion, collision, or long-range 
coulomb interaction, could slightly deplete the Maxwellian tail at high energy has been 
evoked in the past. A detailed study showed that a very small deviation of 0.5 % of this 
Maxwellian distribution leads to a change in reaction rate of 5% in pp reaction, -36% 
in {^He,^ He), -50% in (jBe,p), and -67% in the CNO cycle. The adverse impact on 



the sound speed (Turck-Chieze et al. , 2001a) coupled with the results from the SNO 



experiment, do not support the idea that the Maxwellian distribution of the particle 
velocities is deformed at high energy. 



4-3.6. Is there some mixing in the very central core? When the neutrino puzzle 
surfaced, the possibility was raised for mixing in the solar core to reduce the 
'^Be composition (which is extremely peaked) together with the ^He content and 
consequently to reduce the central temperature and the emitted neutrino flux 



(Bahcall, Bahcall & Ulrich, 1968 Schatzman et al., 1981 Lebreton & Maeder 1987). 



This idea was supported by the fact that the central solar conditions are not far from 
convective instability. Such an effect has an impact on the whole radiative zone far 



greater than what is suggested by the observed acoustic modes (Turck-Chieze et al 



2001a), see figure 3 of Turck-Chieze et al. (2001a). However some small and localised 



mixing in the deep core could be attributed to an initial higher mass (see below). 

A similar idea that the central solar core may be cooled down by the presence of 
WIMPs (Weakly Interacting Massive Particles) acting as a conductive medium, and 



transfering heat from the very center to the rest of the radiative zone (Faulkner, Cough 



and Vahia, 1986), was emitted. More that ten years ago, we indeed used seismology to 



show that such an idea is not favored (Kaplan et al. , 1991 ). At that time the constraints 



on the solar core were not as stringent. The limits imposed by the Sun are regularly 



re-estimated (Lopes, Bertone & Silk, 2002). A revision could be justified by the recent 
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cross-section updates and by the inclusion of the new abundances, but it is probably 
better to wait for the inclusion of the impact of the gravity modes on the density and 
sound speed profiles. The difficulty to separate the effect of mass loss from the effect of 
dark matter hampers progress on this issue (also see section 8). 

In conclusion, it was shown in this section that the physics of the solar interior has 
been greatly improved during the last ten years and that at each stage the predicted 
neutrino fluxes were compared to detection results. The most obvious impact is on the 
neutrinos whose predicted flux changed by more than 50% during that period, as 
shown in Table 6. In parallel, the neutrino flux prediction of Bahcall and collaborators 
varied between 7.5 and 5.8 xl0^cm~^s^^. After the launch of the SoHO spacecraft, 
the predictions of the two teams converged to comparable numbers when using the 
same physics. Since this period, the improvements of the SSM are systematically 
confronted to hehoseismic observations. This second probe was also useful to help 
rejecting several non-standard ideas impacting on the neutrino flux predictions, and 
to qualify the standard solar model which only introduces the main physical processes 
describing the evolution of stars. 



5. The seismic view of the solar interior and comparison with solar models 



In this section, we briefly describe the hehoseismic tools. Other reviews can be read 



for details: Turck-Chieze et al. (1993); Leibacher et al. (1999); Christensen-Dalsgaard 



I (2002); Antia & Basu (2007). Acoustic waves are generated by the granulation at 
the solar surface and propagate inside the Sun to a depth depending on their initial 
velocity. These waves generate very small motions in the solar atmosphere that are 
detectable. We formally treat this information through a perturbation theory because 
the amplitudes of these perturbations are small and the Sun is, to a good approximation, 
spherical. 



5.1. The formalism 

The Sun, as a self-graviting sphere of compressible gas, oscillates around its equilibrium 
state with a period of about 5 min. These oscillations are interpreted as a superposition 
of waves propagating inside the star (acting as a resonant cavity), and forming 
standing waves: the eigenmodes of vibration. By projecting these modes onto spherical 
harmonics YJ™, we write any scalar perturbations as (in the case of the Eulerian pressure 
perturbation p) following Christensen-Dalsgaard & Berthomieu (1991): 

p{r, 9, ip, t) = p'{r)Yi'^{9, (p) exp iujn,i,mt 

and the ^ displacement vector as 

ar,0,V,t) = [^rir),Ur)^,Ur)—^^ l^"^(e, ^) exp(^u;„,,,„t) (5.1) 

where = l/('^^^)[p'/P + is the horizontal displacement, $' the gravitational 
potential perturbation, Un^i^m the eigenfrequency, and p the gas density. The quantum 
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Figure 4. Eigcnfunctions of an acoustic mode {I =0 radial mode, n=23) (left panel), 
and of a gravity mode (£ =2, n=-10) (right panel). 



numbers n, Z, m are respectively the radial order (number of nodes along the radius), 
the degree (the total horizontal wave number at the surface is kh ~ L/Rq, with 
L = \/l{l + 1)) and the azimuthal order (number of nodes along the equator with 
\m\ < I). Restricting the phenomenon to adiabatic oscillations within the Cowling 
approximation ($' neglected) and considering only small radial wavelengths compared 
to Rq (the solar radius), the 4th-order system equations are reduced to second-order 
wave equations, with the following dispersion relation: 




(5.2) 



where the squared norm of the wave vector is written as the sum of a radial and 
a horizontal component |A;| — k"^ -\- A;^, A;^ = l^s horizontal wave number, 

Ff = L'^c^/r'^ is the Lamb frequency, A^^ = g[l/ridlnp/dr — d\np/dr] is the Brunt- 
Vaisala frequency, cu^ = c^(l — 2dHp/dr)/AH'^ is the acoustic cut-off frequency (~ 5.8 
mHz), = —dlnp/dr is the density scale height, Fi is the first adiabatic exponent, 
and = Tip/ p is the sound speed. 

The oscillatory solutions to the wave equation define two types of waves (and 
modes, since some of these waves constructively interfere), namely acoustic ones (with 
^n,i,m > N,Fi) and gravity ones (with ujn,i,m < N,Fi). Figure 4 shows the properties 
of these waves: the acoustic modes have their maximum amplitude at the surface (left 
side), while the gravity modes are excellent probes of the solar core and are evanescent 
at the surface with a rather small amplitude (right side). 

About 3500 acoustic modes (the so-called 5-min oscillations) have already been 
observed. A refined analysis of their properties allows the obtention of a stratified 
information about the solar internal structure from the surface to the solar core through 
the sound speed profile. The gravity waves remain the best probes of the region of 
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Figure 5. Evolution of the sound speed during the hydrogen-burning phase 



neutrino emission and also provide information on the rotation in the core. 

There are two ways to apply the seismic data to probe the internal structure of 
stars: 

a) The direct method: comparison of predicted with observed acoustic wave frequencies; 

b) The indirect method: using inversion procedures to deduce the solar radial profile 
of fundamental variables like the squared sound speed c^, the density p, the adiabatic 
exponent Fi, or the rotation rate, and compare them to a computed model. 

Since the launch of SoHO, we have been using two variables obtained from the 



space instruments GOLF: Global Oscillations at Low Frequency ( Gabriel et al. 1995 ) 



and MDI: Michelson Doppler Imager (Scherrer et al. , 1995), namely the sound speed 
profile and the internal rotation profile, and compared them to up-to-date solar models. 
It allowed an improvement of these models when the agreement was not satisfactory, or 
the introduction in these models of some extra physical phenomena to check whether the 
agreement could be improved. Indeed, a lot of physical phenomena have very specific 
signatures which are informative. We could also potentially access data on the magnetic 
field, but until now we only managed to put upper limits on the field strength due to 
the difficulty to accurately derive such a quantity. 



5.2. The sound speed: a very useful hut demanding quantity 

Since the launch of SoHO in 1995, the solar sound-speed profile has been derived 
with more and more accuracy from helioseismic data. Consequently, if it could be 
demonstrated that there is a significant discrepancy between the observed Sun and 
the SSM, this would be an opportunity to reach beyond the classical theoretical solar 
model. Moreover, in the meantime, this helioseismic information is used to deduce 
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neutrino fluxes constrained by the seismic observation of tfie Sun. 

Before commenting on tlie role of the seismic probe, it is interesting to compare the 
sensitivity of the two probes at hand, the neutrino fluxes and the sound-speed profile, to 
the conditions at the center of the Sun. We reminded the neutrino flux dependence 
on the central temperature in table 6. It is noteworthy that if the Sun was exactly at the 
beginning of the hydrogen-burning phase, the prediction for the neutrino flux detectable 
by the chlorine experiment would be 0.57 SNU (instead of the roughly 7 SNU currently 
estimated) while for the flux detectable by the gallium experiments it would be 67 SNU 
(instead of 127 SNU), meaning a factor, respectively, 15 and 2 smaller. 

On the other hand, as is shown on figure 5, the central sound speed only varied by 
9 % during the same timespan. In fact, as 

AcVc^ = AT/T - A/i//i (5.3) 

and AT/T = 13.5% {Tinit= 13.5 x lO^/T), we have A/i//i = 32 % (the mean molecular 
weight varies from 0.31 to 0.41), and therefore Ac/c= -9%. Consequently the sound- 
speed profile needs to be known with a high relative precision (better than 10~^) in order 
to detect a temperature deviation of less than 1%. This is exactly the precision required 
to test the validity of the present solar structure (see table 7 listing the sensitivity of 
the sound speed to different ingredients of the solar model). An uncertainty of 2 % on 
the central temperature leads to a reduction by a factor of 2 on the neutrino flux. 
This challenge needs to be addressed by the SuperKamiokande and SNO experiments. 
On the helioseismic side, the relative precision of the measurements is very high (about 
10~^) because those are metrological measurements. Therefore helioseismology offers 
very interesting constraints provided that the theoretical acoustic-mode characteristics 
can be validated. The accuracy of the oscillation frequency determination is given by the 
duration of the observation (which can be years or decades) if one can verify that these 
frequencies do not change with time. Their stability was demonstrated for frequencies 



below 1.6 mHz (Couvidat, Turck-Chieze & Kosovichev 2003 Garcia et al. 2004) with 



the GOLF experiment. 

Contrary to boron neutrinos, acoustic modes do not provide a direct determination 
of the temperature. Ac^/c^ does not depend only on AT/T as was previously shown. 
Nevertheless, if we have strong constraints on the density and pressure, we also have 
an indirect but strong constraint on the central temperature. In 2001, we investigated 
all the different issues mentioned in table 6 and we reached the required precision to 
predict neutrino fluxes constrained by helioseismology. These fluxes were compared to 



the very important SNO results (Turck-Chieze et al. 2001b, 2004a Couvidat, Turck- 



Chieze & Kosovichev, 2003). During these years, major efforts have been carried out 
to improve the quality of the seismic indicators because the acoustic mode frequencies 
are largely influenced by the physics of the outer solar layers. The different techniques 
used contributed to identify biases in the oscillation frequency determination. Three 
phenomena have been thoroughly studied for investigating the physics of the solar core: 
the stochastic excitation of acoustic waves, the influence of the solar activity on the 
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Table 7. Sensitivity of the sound s 


peed to the physical processes 


Quantity 


variation 


Ac^/c^ variation 


T 


1 % 


1% 




1 % 


0.1 % 




4 % 


0.1 % 


X^He 


25% 


0.1 % 


(p,p) reaction rate 


1% 


± 0.1% 


(^He, ^He) reaction 


- 25 % 


- 0.1 % 


rate 






(^He, ^He) reaction 


-25% 


+0.2% 


rate 






(p, ^Be) reaction rate 


10% 


none 


(p, ^^0) reaction rate 


-50% 


- 0.1-0.2 % just at the center 




Figure 6. left: Relative difference between the squared sound speed in the Sun 
extracted from GOLF+MDI and the squared sound speed of reference solar model in 
function of radius, compared with a model which does not contain the microscopic 
diffusion (dashed line). From Brun, Turck-Chieze & Morel (1998). 



absolute value of the mode frequencies, and the asymmetry of the mode distribution 



due to the interaction of the modes with the solar background (Basu et al. , 2000 



Couvidat, Turck-Chieze & Kosovichev, 2003). SoHO/GOLF detected modes of low 



frequency which have a long lifetime and propagate inside solar cavities for which the 
noise coming from stochastic excitation of waves and the effects of the solar cycle are 
both reduced. The relative accuracy of the sound-speed determination in the central 
solar region was improved by an order of magnitude and is now of the order of 10~^ in 
the core (Turck-Chieze et al. , 2001b Couvidat, Turck-Chieze & Kosovichev, 2003). 
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Figure 7. Comparison of the sound speed in the core deduced from the acoustic 
modes and the sound speed of standard models and seismic models (in black). 
Acoustic modes explore the solar interior down to 6% so reasonably well in the 
region of emission of boron neutrinos. Left: standard model of the nineties, 
Right: new standard model, (points with error bars joined by straight lines). 
Adapted from Turck-Chieze, Piau & Couvidat (2011). 



5.3. The sound-speed and density profiles: observations and predictions 



Figure 6 shows the first squared sound-speed difference obtained using GOLF and MDI 
data aboard SOHO, between a standard model including or not the gravitational settling 
of elements. Figure 7 shows a zoom on the solar core of the sound-speed profile of 
standard models (obtained with old and new abundances) and of the seismic model 
(see below). The vertical error bars are so small that they are almost invisible on the 
figure: such a high precision is due to the long temporal duration of the observation 
considered, and to the metrologic character of the Doppler velocity method (assuming 
that there is a good control of the instrumental noise). The horizontal error bars 
translate the sensitivity of each mode to the physics of the solar core. The large number 
of points compensate for these large error bars. The values of the sound speed in the 



core have been recently confirmed with ground-based observations (Basu et al. , 2009) 



The development of helio seismology (also see section 1) and the recent progress on CNO 
photospheric abundances, lead to the conclusion that the SSM reached its limits when it 
comes to representing the real Sun. 

Figure 8 shows the evolution of the discrepancy between the observed sound speed 



and density and the SSM predictions. Today, a difference between SSM, for Caffau et 



al. (2009) and Asplund et al. (2009) compositions, and seismic observations still exists. 



which is much larger than the seismic error bars. 



5.4. The gravity modes 

Gravity modes are the best probes of the solar core. They inform us on the density and 
rotation profiles of the core and consequently help us connect the current state of the 
central solar core to its past states and to the physical quantities characterizing them. 
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Figure 8. Relative differences between (a) tlie square of the sound speed and 
(b) the density deduced for the Sun using the GOLF/MDI frequencies and 
those of different Saclay models. Left: Updated standard models including 
new estimates of CNO abundances by Asplund et al. ( 2004[ ) (full line with 
error bars coming from seismic observations: tac A model, dot line: tac L 
model, dot dashed line: tac H model) and the seismic model (full line). From 



Turck-Chieze et al. (2004a). Right: New SSM using Asplund et al. (2009) 



composition and CGM convective prescription. From Turck-Chieze, Piau & 



Couvidat (2011) 



Gravity modes are been searched for by the three hehoseismic instruments onboard the 
SoHO satelhte. The Doppler velocity technique, with measurements performed in the 
solar atmosphere where the noise is smaller than at the surface, proved the most adapted 
to this search. Indeed, the amplitudes of gravity modes are small and these modes are 
in the frequency range 10 to 400 /iHz, which is dominated by granulation noise. In 
fact, two years after the launch of SoHO, the GOLF instrument, which looks at the Sun 
globally, appeared to be the most promising instrument for this search, thanks to its 
very low instrumental noise. Two publications presented positive detections of patterns 
in the gravity mode range: the first publication was dedicated to the range of mixed 
modes (modes that present characteristics of both gravity and acoustic waves) where 



two patterns were followed over a period of 10 years (Turck-Chieze et al. , 2004b Garcia 



et al. , 2008 ). These patterns could favor an increase in the rotation rate in the solar core. 
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but no clear determination of the central frequency of a potential quadrupole mode (n= 
-3) was obtained so far. The second paper is dedicated to the frequency domain where 
gravity modes are equally spaced in period. By cumulating the power of 20 oscillation 
modes, we highlighted a pattern attributed to dipole modes (Garcia et al. , 2007). Some 
of these modes are now analyzed individually and at least 6 modes are properly identified 
(in the range n = -4 to n= -10). These results should be published this year. A definitive 
detection of these gravity modes improves the sound speed, density and rotation profiles 
in the core. The rotation rate of the very central region already seems larger, by a factor 
5 to 7, than the rest of the radiative zone (figure 11) which is a very important result 
of the SoHO spacecraft. The frequencies of the detected modes are very close to the 
values predicted by the seismic model in Couvidat, Turck-Chieze & Kosovichev (2003); 



Mathur et al. (2007). 



6. The detection of solar neutrinos and comparison with predictions 

Solar neutrino detectors focus on neutrino energy lower than 20 MeV, through one or 
more of the three different types of reactions : 

- neutrino capture on a nucleus (A,Z) by W Exchange Process (WEP)[||j 

+ (A, Z) ^ X- + (A, Z + 1) (6.1) 

Experiments based on this reaction are only sensitive to z/e . When v^. oscillate and 
become or Vr , they cannot be detected anymore by this process since the threshold for 
producing a muon (m= 106 MeV) or a tau (m= 1784 MeV) is well above the maximum 
solar neutrino energy (~ 14 MeV). An electron is (almost) isotropically released by the 
reaction involving z/g , thus giving (almost) no information on the neutrino incidence 
direction. However the electron energy spectrum directly reflects the neutrino energy 
spectrum: Eg = Ej^ - Ethreshoid where Ethreshoid is the energy threshold of the reaction. 

- neutrino scattering by Z exchange process (ZEP): 

i/x + A ^ z/, + A* (6.2) 

The reaction cross section does not depend on the neutrino flavor. This type of reaction 
is very useful to distinguish neutrino oscillations (constant flux but mixed flavors) from 
a genuine neutrino flux deficiency. 

- elastic scattering on electrons: 

z/x + e" -)■ z^x + e~ (6.3) 

This reaction occurs both via WEP and ZEP for z/g , and only via ZEP for z/^ and v-r . 
The resulting cross section for the scattering is about 6 times larger than for the z/^ or 

II The denominations "charged current" process (CC) and "neutral current" process (NC) have been 
replaced here by "W exchange" process (WEP) and "Z exchange" process (ZEP) which better describe 
the physical mechanisms. 
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Ur scattering. This means that the elastic scattering of Ve on electrons proceeds mainly 
via WEP. An advantage of this reaction is that, for kinematic reasons, the scattered 
electron direction is correlated with the direction of the neutrino. This property is 
very helpful for background noise reduction. However the electron energy spectrum is 
not identical to the neutrino spectrum, though scattering experiments give significant 
information on this neutrino energy spectrum. 

A summary of the solar neutrino experiments is given in table 8. 



6.1. Radiochemical experiments 

These experiments use the method of neutrino capture on a nucleus (A,Z) by W exchange 
process (WEP) and detect only . 



- Homestake experiment: 1968-1995 

This historical experiment follows the initial idea of Pontecorvo (1947). It considers 
the reaction z/e + '^^Cl — )■ ^^Ar + e~ to capture solar ■ The energy threshold is 
0.814 MeV, and the resulting radioactive ^^Ar isotopes decay by electron capture with 
a half-life of 35 days. The cross section varies from 5x 10~^^ cm^ at 1 MeV to 10~^^ cm^ 
at 14 MeV. Therefore, about 10^^ atoms of chlorine are necessary to observe one 
atom a day. Consequently, from several hundred of tons of a compound of chlorine, one 
atom of ^''Ar a day (or a few atoms every other month) needs to be extracted, and its 
decay needs to be unambiguously observed. 

A big tank containing 615 tons of perchlorethylene C2Cl4was settled in 1967 in 
the Homestake gold mine (South Dakota, USA), at a depth of 4100 ±200 hg/cm^of 
standard rock to shield against cosmic rays. The first results gave an upper limit of 
3 SNU ( Davis, Harmer & Hoffman , 1968 ) , more than a factor 2 lower than the predictions 
from solar models. 

A run consists of three steps: the exposure, the argon extraction, and the counting 
of the ^^Ar atoms. Before each run a small amount of carrier gas (about 1 cm^of non 
radioactive argon isotopes, ^^Ar or '^^Ar) is stirred and dissolved into the C2CI4. The 
recovery of this carrier allows a measurement of the argon extraction efficiency. The 
exposure to solar neutrinos lasts about two months. Argon (the carrier and the ^'^Ar 
produced by solar z/e or any background), is then removed from the C2Cl4by circulating 
a few hundred cubic meters of helium through the tank. A charcoal trap cooled at 
liquid-nitrogen temperature absorbs the argon and lets the helium pass through. The 
argon is then desorbed and pushed into a small proportional counter (about Icm^). 
10 % of methane is added to optimize the working of the counter. The counters are 
made of low-activity materials and put in a lead shielding located in the mine itself. 
^''Ar decays by electron capture with a half-life of 35 days and emits 2.82 keV Auger 

^ Underground detectors are located at different depths of different kinds of rocks. In this field, one 
generally gives the depth with a unit called hg/cm^of standard rock. Standard rock is defined as rock 
with a density of 2.5 g/cm'^, a mean atomic number Z—11, and a mean atomic mass A— 22. 
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detection 


signature 




experiment 


reaction 


threshold 


main contributors 


events /yr 


Homestake 


z/e + 3^C1 ^ 37Ar+ e- 


0.814 MeV 


37Ar 


100 


615 tons 










1968-2003 










Kamiokande 


z^e e~ — >-i/e e~ 


5 MeV 


E > 7.5 MeV 




2140 tons H2O 










1987-1995 






Sun direction 




SAGE (1990- 


i/e + ^^Ga ^ ^^Ge + 


0.233 MeV 


30 ^^Ge 




50 tons 










Gallex (1991-1997) 






for the 3 detectors 


20 


30 tons 










GNO (1998-2003) 






t'pp ; t^Be , ^ , 


20 


radiochemical 










Super-Kamiokande 


z^e e~ — >-i/e e~ 


5 MeV 


E > 5 MeV 


24000 


50000 tons H2O 










Cerenkov 






Sun direction 




1999- 










Sudbury 


Ve D — >e~ p p 


6.5 MeV 


E > 6.5 MeV vb 


9750 


1000 tons 






(D2O) - (H2O) 




D2O 


i/x D -^Uy^ p n 


2.2 MeV 


n capture on "^^Cl ub 


2800 


Cerenkov (1996- 






(D2O + NaCl) - (H2O) 




(2000- 


z^e e~ — >i/e e~ 


5 MeV 


E > 5 MeV z/B 
Sun direction 


1100 


Kamland 


z^e + P — ^e"^ -|- n 


1.8 MeV 


reactors 




1000 t liquid scintillator 










(2002- 










Borexino 


z/e e~ — >-i/e e~ 


0.2 MeV 


E > 2.8 MeV ub 


130 


2000 tons scintillator 






0.2< E < 0.8 MeV z/Be 


11000 


(278 t ) 






no Sun direction 




(2007- 


z^e "B ^"C* e- 




E > 3.5 MeV 
no 7 or 7(2,4.3,4.8 MeV) 


2300 




z^x ^^B ^^^B* i/x 


4.5 MeV 


7 










(4.4 or 5 MeV) 




Icarus I 


z/e ^°Ar ^40k* c- 


11 MeV 


E > 5 McV 


300 


600 tons 






7 2.1 MeV 




liquid argon 


z/x ^"Ar ^4"Ar* z^x 


6 MeV 


7 


60 


drift chamber 






6.1,7.8,9.6 MeV 




May 2010 


z/e e~ -^Ve e~ 


5 MeV 


E > 5 MeV 
Sun direction 


240 
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electrons. The counter background noise is reduced by a pulse-shape analysis of the 
signal: ^^Ar decays show a characteristic short rising time. The counter background 
noise has been reduced to 0.01 count/day after several years. The average ^^Ar count has 
been of about 0.5 atoms/day. The major source of background noise comes from cosmic 
ray muons (about 5 /m^/day at this depth). They produce low energy protons either 
directly or through the hadronic cascade which may follow. These protons produce ^^Ar 
through (p,n) reactions with ^^Cl . Two indirect determinations of this background lead 
to a value of (0.08 ±0.03) '^'^Ar /day (about 20% of the signal as will be seen below). 
The other sources of background noise, fast neutrons from surrounding rock and cosmic 
ray neutrinos, have been found negligible. 

More than one hundred runs have been performed since 1968, almost continuously 



except in 1985-1986 (Davis, 1994). The first two years of runs are no longer considered, 
due to considerable improvements and modifications during this period. The final result 
for 25 years corresponds to a detection of 2200 atoms of '^^Ar. The average solar neutrino 
flux of the 108 runs covering the period 1970-1995 is 2.56±0.16(stat.) ±0.16(syst.) SNU 



(Cleveland et al. , 1998). R. Davis received a joint Nobel Prize for this work, with M. 



Koshiba for the Japanese research on this subject (see below). 

- Gallium experiments 

The main objective of the radiochemical gallium experiments is the detection of 
the z/pp which are produced in the primary pp fusion reaction using the fact that the 
reaction z/e +^^Ga — J-^^Ge +e~ has an energy threshold of 233 keV, significantly below 



the maximum value for i^pp (420 keV). This is an original idea of Kuzmin (1966). The 
^^Ge decays (half-life of 11.43 days) by electron capture (the probabilities of K, L, and M 
captures are respectively 87.7%, 10.3% and 2%), emitting low-energy Auger electrons 



and X-rays. Two gallium experiments have been designed: SAGE in Russia (Gavrin 



et al. 1990) and GALLEX in Italy (Anselmann et al. 1992). Their principle is the 



same: the gallium is exposed to solar neutrinos in a low-background environment, then 
the ^^Ge atoms produced by the reaction are extracted by a chemical method and 
transformed into a counting gas (germane GeH4 ). A proportional counter counts them. 
The main difference between these two experiments is that GALLEX uses a solution of 
GaCls as the target while SAGE directly uses the gallium metal. 

- Gallex and GNO: 1991-2003 

The Gallex experiment, a Germany-Italy- France USA-Israel collaboration, was 
installed in the Gran Sasso Underground Laboratory (Italy), in a highway tunnel which 
crosses the Apennine mountains. The detector was a cylindrical tank (8 m high and 
3.8 m diameter) containing 30.3 tons of gallium in the form of GaCls (8.2moles/liter of 
GaCla and 1.9 moles/liter of HCl). The expected number of ^^Ge atoms was 1.2 per day 
according to theoretical predictions, but in fact it has been less than 1 / day. 

Every three weeks a large flow of nitrogen (150m^/h) circulated through the tank 
for about 20 hours to remove the germanium. Germanium, in the form of GeCl4 , is 
very volatile in presence of HCl and was taken out with the nitrogen. It was absorbed 
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into pure water inside a system of exchange columns full of glass helices which ensured a 
large surface of contact. A concentration step was then followed by a transformation of 
the GeCU into germane GeH4 . A small amount of a given stable isotope of germanium 
(''^Ge, "^^Ge, '''^Ge, "^^Ge), was added as a carrier in the tank before each run. It allowed 
a measurement of the extraction efficiency and was used as the counting gas. The gas 
mixture in the proportional counters was germane (30%) and xenon (70%). The small 
proportional counters (Icm^) were made with ultrapure material. Only K (10.4 keV) 
and L (1.25 keV) electron-capture peaks could be observed, with an energy resolution 
of about 20 %. The first signature of a ''^Ge decay was a count observed in the K or L 
peaks. Backgrounds as low as 0.06 count per day (respectively 0.01) have been obtained 
for the L-peak (respectively K-peak). The extraction efficiency was larger than 95% 
and the counting efficiency was estimated for each counter (about 65%). In a 21-day 
run, GALLEX observed about 3-4 ^^Ge decays taking into account all the efficiencies 
and natural decays. 

The main source of background noise was the reaction ^^Ga (p,n) ^^Ge . The 
protons came either from natural radioactivity of the liquid solution or from cosmic ray 
muon flux (about 15 /m^/day at a depth of about 3500hg/cm^) or from neutrons from 
the surrounding environment (building, rock, tank, etc.). The purity speciflcations 
required for the solution have been fully satisfled (0.03 pCi/kg of ^^^Ra, less than 0.05 ppb 
of U and less than 0.1 ppb of Th) and induced a contamination < 0.1 SNU. The protons 
coming from the interactions of the cosmic ray muons induced a signal of about 3 SNU. 
The fast neutrons in the environment (> 1 MeV) reacted through n -|- ^*^Ca — )■ ^^Ar -|- a. 
The induced background was estimated to be about 0.5 SNU (Cribier et al. , 1995 ). The 
total background was < 4 SNU. Solar neutrinos only produce ^^Ge from ^^Ga but 
cosmic rays produce ^^Ge (half-life of only 1.63 days) and also ®^Ge which decays to 
^^Ga by electron capture with a half-life of 288 days (^®Ga , whose lifetime is 68 minutes, 
decays mainly by positron emission). Thus, its decay is identical to the ^^Ge decay. The 
gallium target was exposed to cosmic rays at the ground level for several months but 
this problem was solved by heating the solution. 

The results of the different runs were presented for different periods of time, Gallex 



1: Anselmann et al. (1992); Gallex 2, 1998: Hampel et al. (1998); GNO 2000: Altmann 



et al. (2000). The flnal value of the joint 123 runs was 69.3 ± 5.5 SNU, confirming the 



deficit of neutrino detection (Altmann et al. , 2005), see also table 9. 
- SAGE 1989- 

The Soviet- American Gallium Experiment is located in the underground laboratory 
(Baksan Neutrino Observatory) of the Baksan Valley in the Caucasus mountains (ex- 
USSR), under about 4700hg/cm^. The muon flux is about 2 /m^/day, i.e. 7 times less 
than in the Gran Sasso laboratory. SAGE started to take data in December 1989 with 
a 30-ton gallium target (Gavrin et al. , 1990). At that time, the counter calibration was 
done with an ^^Fe source. The data analysis is based on a maximum likelihood method 
including ^^Ge decay and a constant background. Identifying no ''^Ge decay (after about 
two months of counting each run), the flrst result of the SAGE collaboration was an 
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upper limit on the solar neutrino flux but with a strong suspicion of deflciency. 

The gallium is in the form of gallium metal (57 tons at the beginning and then 49 
tons since 1998) placed in reactors of 7 tons each, equipped with stirrers and heaters 
that maintain the temperature just above the melting point (29.8 C). The advantage of 
using the metal is its reduced sensitivity to both internal and external backgrounds and 
a smaller volume. The major disadvantage is the chemical treatment. The chemistry for 
the extraction procedure is a little more complicated than for a liquid GaCla solution, 
(Gavrin et al. , 1992) for details, and requires additional steps and somewhat greater 
quantities of fresh reagents for the extraction. The efficiency is well known and of 
the order of 80-90% instead of the 95% reached in the Gallex experiment. Low level 
radioactivity materials have also been used in this experiment. The cosmic muon 
induced background is negligible, being < 0.5 SNU. 



A test of the efficiency of the overall procedure has been performed (Abdurashitov 



et al. , 1996). An artificial neutrino source (using 200 g of highly enriched chromium 
: 86% of ^°Cr ) of radioactivity > IMCi of ^^Cr , was used to calibrate the detector 
and to obtain an overall consistency check of the whole experiment. ^^Cr is obtained 
by irradiating ^°Cr with thermal neutrons. The source was produced from enriched 
chromium (40% of ^"^Cr instead of 4.5% in the natural chromium) in the Siloe reactor 
located in Grenoble, France. 

Assuming the solar neutrino production rate was constant during the period of 
data collection, the combination of 168 extractions from December 2007 until the end 
of 2008 returns a capture rate of solar neutrinos with energy higher than 233 keV of 
65.4 — 3.0 + ?>.l{stat) — 2.8 + 2.Q{syst) SNU. The weighted average of the results of 
all three Ga solar neutrino experiments, SAGE, Gallex, and GNO, is now 66.1 ± 3.1 
SNU, where statistical and systematic uncertainties have been combined in quadrature. 
During a recent period of data collection a new test of SAGE was made with a reactor- 
produced ^''Ar neutrino source. The ratio of observed-to-calculated rates combined 
with the measured rates in the three prior ^^Cr neutrino-source experiments with Ga, 
is 0.87 ± 0.05. A proton-proton flux of 6.0 ± 0.8 10^°/cm^/s was derived. 



6.2. The real time experiments: Kamiokande, SuperKamiokande, SNO a,d KAMLAND 

These experiments have contributed to signiflcantly increase the statistics of the solar 
neutrino detections. The flrst experiment, Kamiokande and, later, Superkamiokande, 
detect the elastic interaction on electrons. The next generation of experiments, Sudbury 
Neutrino Observatory, Borexino, and Icarus, aim to measure the contributions of both 
the W and Z exchange processes and to detect the elastic interaction on electrons. 
- KamiokaNDE, SuperkamiokaNDE experiments: 1986-1995 , 1999- , 2002- 

The flrst real time experiment was the Kamiokande experiment. Its principle was 
based on elastic neutrino scattering: i/e + e~ — )■ z/g -|- e~. The scattered electron is 
detected through the Cerenkov light emitted and its direction is strongly correlated 
with the direction of the incoming neutrino. The target consisted of a large cylindrical 
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tank (16 m high and 15.6 m diameter) containing 2140 tons of ultra-pure water. The 
Cerenkov hght was captured by about 1000 photomultiphers of 20 inch diameter covering 
20% of the wall surface. This detector, located in a mine at a depth of 2700hg/cm^of 
standard rock (Hirata et al. , 1988), was primarily designed for observing the nucleon 
decay with a lower limit of 2.6 x 10^^ years. It was then surrounded by a water Cerenkov 
anticoincidence layer of thickness 1.4 m and started to collect solar neutrino data at the 
end of 1986. Its great success was the observation of 12 neutrinos from the supernova 
SN1987A on February 23, 1987. 

The main background sources were the radioactivity in the water (/3 — 
decay of ^^^Bi) (E< 9MeV), 7 rays from surrounding rock, and unstable spallation 
products by cosmic muons (E> 9MeV). These backgrounds do not have any correlation 
with the Sun's direction. The water Cerenkov anticoincidence layer absorbed 7 — rays 
from surrounding rocks and monitored the cosmic muons. The useful volume (called 
fiducial volume) for the solar neutrino detection was reduced to 680 tons (about 2.5 m 
from the wall of the detector), to enhance the signal-to-background ratio. The electron 
detection efficiency was 50 % (respectively 90 %) at 7.6 MeV (respectively 10 MeV) in 
1987 and has continuously improved since. The selected event sample needed to have a 
directional correlation with the Sun. 

The enhancement around cos 9 = 1 (where 9 is the angle between the trajectory of 
the electron and the direction of the Sun) corresponded to solar neutrinos. It constituted 
a direct evidence for a detection of solar neutrinos, but the neutrino energy spectrum 
and the experimental electron spectrum both needed to be known to determine the 
neutrino flux. The theoretical prediction (including Monte-Carlo simulation and the 
energy spectrum calculation) was deduced from Bahcall & Ulrich (1988). The first 
result of Kamiokande indicated that the observed signal was less than expected from 
solar models. 

After nine months used to improve the detector, new data were taken from January 
1991 onwards. The combined value accumulated from 1986-1991 is 0.49 ± 0.05 ± 0.06 



dTotsuka et al.[|199lD times the prediction of [Bahcall k Ulrich] ( |1988| ) and 0.74 ±0.07 ± 
0.09 times the prediction of Turck-Chieze et al. (1988). The trigger rate was about 0.5 
per second and the expected solar neutrino signal about 0.5 per day above 7.5 MeV! 
Consequently, it was not possible to assign an individual event to a solar neutrino 
interaction. 

SuperKamiokande is a significant extension of Kamiokande. It uses a tank of 50,000 
tons of pure water providing 65 solar neutrinos per day in a 22,000 tons fiducial volume 
(Hirata et al. , 1988). The SuperKamiokande experiment was ready in April 1996. 
Searches for possible day-night and semi-annual variations of the ^B solar neutrino 
flux have been carried out, using the 1040 days of Kamiokande II data. No such short 



time variations were observed (Hirata et al. , 1991). Taking into account the various 



uncertainties, this results in an upper limit of about 30% on the relative difference in 
the day/night neutrino fluxes. The detection results obtained with SuperKamiokande 
contributed to the Nobel Prize of 2002, partly awarded for the demonstration of the 
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existence of neutrino oscillations. 

After full restoration of the detector in 2006 (following implosions of 6,600 
photomultipliers in 2001), the energy threshold has been lowered down to 4.5 MeV. 
The mean values of the Superkamiokande measurements for all three experiment phases 



are respectively: 2.35 ± 0.08,2.38 ± 0.17,2.31 ± 0.0510Vcm7s (Hosaka et al. 2006 



Cravens et al. , 2008). The day-night effect, which is a direct test of matter effect, does 
not clearly show any asymmetry within a 2-3% error bar, and this error should be further 
reduced to 1.5%. 

- The Sudbury Neutrino Observatory experiment: 2000-2006 

The Sudbury Neutrino Observatory (SNO) experiment (a Canada-USA-UK 



collaboration) (Aardsma et al. , 1987) consists of 1000 tonnes of heavy water 
D2O surrounded by 5 m of purified light water H2O . The Cerenkov light emitted 
by electrons is detected by photomultipliers, like in the Kamiokande experiment. The 
project was constructed in a deep mine (2070 m underground) near Sudbury, in Canada, 
and began taking data in 2000. The main advantage of this experiment is its ability to 
detect both z/e (through the W exchange process on deuterium) and all neutrino types 
(through a Z exchange process producing dissociation of deuterium into a neutron and 
a proton). The WEP reaction is detectable via Cerenkov light from electrons which are 
produced with an energy nearly equal to the incoming , thereby providing about 20 % 
energy resolution to search for possible distortion (through the MSW effect) of the 
energy spectrum. The ZEP reaction is detected via the 8 MeV 7 rays produced when 
the neutron is captured in ^^Cl from NaCl cycled in and out of the heavy water. Discrete 
^He proportional counters have been developed as an alternative neutron detection 
technique. 

The major experimental challenge is to reduce the backgrounds to a very low level. 
This requires the use of low activity materials: less than 10~^^ g/g of U and Th in the 
heavy water; the high energy 7 rays from the U and Th chains can photodissociate the 
deuterium, emitting a neutron which can simulate a Z exchange process. 

This experiment is sensitive only to solar u-q coming from the decay, but presents 
a lot of advantages: measurement of the neutrino spectrum, sensitivity to ZEP reactions, 
and sensitivity to a day-night variation of the solar flux. SNO was key to proving the 
existence of neutrino oscillations by estimating the three neutrino fluxes (electron, muon, 
and tau neutrinos). SNO was measuring both neutrino-electron elastic scattering (ES) 
like Superkamiokande (that means z/e -|- 0.16 Ur -1-0.16 Ur) and the charged-current 
interaction (CC) on deuterons which is sensitive to Ue. Therefore it measured for the 
fist time the sum of all three species, which is directly equal to the number of emitted 
by the Sun. Consequently a direct comparison between predictions of SSM or seismic 
solar model (SeSM) and detection results is possible. In June 2001, SNO obtained a 
detection of 4.95 ±0.72 lO^cm"^ for a prediction of the SeSM of 5.44 ±0.99, which was a 
wonderful success for the two disciplines of neutrino physics and helioseismology. Such 
result triggered the attribution of the Nobel Prize to the pionnier experiments done in 
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Japan and the USA. It also shed a new hght on the SK resuhs which include 0.16% 
of the two other neutrino species (tau and muon), and explained the strong electronic 
neutrino deficit detected by the chlorine experiment. 

SNO measured the neutral-current (NC) disintegration of deuterium which is 
equally sensitive to all active neutrino favors: 

z/x + D-)-z/x + n + p 

SNO captured the neutron produced by the NC reaction. In phase I, the neutron was 
captured mainly by the D2O of the solution. In phase II, the addition of 2 tonnes of 
NaCl to the 1000 tonnes of heavy water enhanced the efficiency of the neutron detection 
by at least a factor two. Moreover the threshold of electron detection for the boron has 
been reduced to 3.5 MeV. In phase III, the neutron detection is also done with an ^He 
proportional counter through the reaction ^He + n — )■ p+^T (983 events per year instead 



of 267 for photomultiplier tubes). The final results published by Aharmim et al. (2008) 
for the SNO collaboration, after the end of the experiment, are : 

•fee = 1.67l°:°l(stat)t°:°3^; $es = 1.77i°:^t(stat)^°:?^; <fNC = 5.54l°:i?(stat)l°:i lOVcmVs 

The charged current is only sensitive to electronic neutrinos, the elastic scattering is 
sensitive to all neutrinos but is dominated by the cross section of the electronic neutrinos 
which is 6 times greater than the cross section of z/^ and z/^, and the neutral current allows 
a count of all the neutrino species. These values are compared to predictions in table 6. 
Altogether, these seminal results demonstrating the detection of all the neutrino species 
and contributing to estimate the oscillation impact on the boron neutrinos, fostered a 
lot of prizes. They also put strong constraints on the neutrino oscillation parameters 
because SNO delivered different fluxes (electronic neutrinos alone or the sum of the 
three flavors, see section 6.5). Finally the electronic-neutrinos boron flux corresponding 



to the total period of detection of SNO is 5.046 lo:i52(stoi)lo:i23 W /cm^/s Aharmim 



et al. (2010b). SNO was turned off on 28 November 2006 even though it also has the 
potential to detect atmospheric and supernova neutrinos, but with a lower rate than 
Super Kamiokande. SNOLAB operates other experiments and prepares the SNO+ 
detector. 



- The KAMLAND detector: 2002- KAMLAND is the largest (1000 t) liquid scintillator 
detector. The objective of this detector is to measure not only the solar neutrino low 
energy spectrum but also to detect anti-neutrinos. It is located in the Kamiokande 
underground mine providing a cosmic-ray muon attenuation by a factor of 10^ with 
respect to the earth surface level. The central spherical vessel of 13 m diameter 
containing the liquid scintillator, is suspended in a transparent buffer oil. The 
photomultipliers (1325+554) cover 34% of vr . An outside tank, flUed with 3200 tonnes of 
purifled water, serves as a radiation shield against neutrons and 7 rays from surrounding 
rocks. Since 2002, the beginning of the data collection, KAMLAND has used nuclear 
reactors to perform neutrino oscillation measurements and geoneutrino detection. The 
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Figure 9. Top: Results of tfef^rsl year of BOREXINO that shows the largest 
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survival probabilities for the different sources of neutrinos and different predictions. 



From Escrihuela et al. (2009) 



first two years have led to 258z/e candidate events with energies above 3.5 MeV, leading 



to a deficit of detection (Haraki et al. , 2005). The distortion of the spectrum puts very 



strong constraints on the oscillation parameters (Figure 10). See section 6.5. 
6.3. Toward a neutrino spectrum: the BOREXINO experiment: 2007- 



The Borex project (Raghavan & Pakvasa, 1988) is another example of real time 



experiment. Nominally, it involves a large tank containing 2000 tons of borated liquid 
scintillator (200 tons of ^^B ) immersed in pure water. It is supported by a USA- 
Europe collaboration in the Gran Sasso underground laboratory (about 3300 m of 
water equivalent) in Italy. A major difficulty is to obtain a very pure liquid scintillator. 
The WEP reaction is marked by an electron in coincidence with a photon. A first 
step of Borex, called Borexino ("il piccolo Borex"), is currently running in the Gran 
Sasso laboratory with a detector of 278 tons of pseudocumene doped with 1.5g/l of a 
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fluorescent dye of dipfienyloxazole (Alimonti et al. , 2009). 



Borexino aims to clarify the nature of the solar electronic neutrino oscillation by 
looking at the ^Be and ®B neutrinos with the same detector. The first objective of 
Borexino is to properly detect the ^Be neutrinos because it is the first experiment to 
report a real-time observation of the low energy spectrum (between 0.2 to 4.5 MeV, see 
figure 9 top left). The nominal counting rate is of 1500 counts per day dominated by 
the residual muon flux. These muons are rejected at the 99% level by an outer detector. 
^''^Tl and ^^^Bi are properly rejected. The pollution by neutrons is reduced to 10~^ 
counts per day. The 0.862 MeV ^Be neutrinos signature is given by the Compton-like 
edge of the recoil electrons at 665 keV. After 192 days, the '^Be flux is estimated at 
49± 3 counts per day per 100 tons with a proper determination of all the contaminants 
(85 Kr, 2iOBi and ^^C). This corresponds to a neutrino flux of 3.36 ±2 ±0.27 x lOVcmVs 



(Arpesella et al. , 2008). Compared to the solar model predictions (see for example 
table 9), the neutrino-fiux deficit is of only about 29%, much less than the 50% (se 
the same table) obtained for the boron flux. This confirmed the energy dependence 
of the survival probability of the electronic neutrinos (see below). Moreover, Borexino 
also detects about 0.26 counts per day per 100 tons of boron neutrinos between 2.8 and 
16.3 MeV, with 0.14 counts between 5 and 16.3 MeV (figure 9 top right). It is the first 
experiment which detects the boron neutrinos with a liquid scintillator detector and a 
threshold of 2.8 Mev, obtaining 75 ± 13 counts. The first analysis gives an elastic flux 



for boron of 2.4 ± 0.4 ± 0.1 (Bellini et al. , 2008, 2010), in good agreement with the 



previous measurements. Borexino also detects the other neutrino sources pp, pep, CNO 
neutrinos and ^^C and ^^B. 



6.4- Neutrino predictions from the standard and seismic models 

The excellent quality of the sound-speed data derived from GOLF+MDI allowed us to 
produce a Seismic Solar model (SSeM) that takes into account this sound-speed profile 
to avoid any possible theoretical bias in the neutrino predictions. Such biases result from 
the hypotheses applied to build the standard model. The seismic model was obtained 
by modifying some physical inputs (mainly reaction rates and opacity) within their 
error bars. The seismic model contributes to stabilizing both the neutrino fluxes and 
gravity-mode frequency predictions. Until a complete Solar Dynamical Model (SDM) 
can be built, the solar seismic model is the best way to deal with dynamical effects that 
are unaccounted for in ID. This model is shown in Figures 7 and 8 and is defined in 



Turck-Chieze et al. (2001b) and Couvidat, Turck-Chieze & Kosovichev (2003). Table 



6 recalls how the boron neutrino prediction has evolved in parallel with the improved 
description of the microscopic physics of the solar models. Table 9 compares the seismic 
model predictions to the results of all the detectors. It is noticeable that the most recent 



prediction of the SSM is not that close to the detected flux (Turck-Chieze et al. 2004a 



Bahcall & Pena-Garay 2004). 



As shown in Tables 6 and 9, the seismic solar model produces flux predictions very 
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Figure 10. Neutrino oscillation contour. Left SNO only; middle: all the solar 
neutrino detectors: SNO, SK, CI, Ga, Borexino; right: Solar neutrinos + 



KAMLAND. From Aharmim et al. (2008). 



Table 9. Neutrino predictions of the present standard solar model and of the seismic 
solar model compared to radiochemical detections and real time experiments (including 
BOREXINO) given in bold characters. See text for the reference of the dilferent 
experiments. 





Predictions without 


Predictions with 




neutrino oscillation 


neutrino oscillation 


HOMESTAKE 




2.56 ± 0.23 SNU 


Standard model 2009 


6.315 SNU 


2.24 SNU 


Seismic model 


7.67 ± 1.1 SNU 


2.76±0.4 SNU 



GALLIUM detectors GALLEX 

GNO 

GALLEX + GNO 
SAGE 

GALLEX+GNO+SAGE 

Standard model 2009 120.9 SNU 
Seismic model 123.4 ± 8.2 SNU 



73.4 ± 7.2 SNU 

62.9 ± 5.4 ± 2.5 SNU 

67.6 ± 3.2 SNU 

65.4 ± 3.3 ± 2.7 SNU 

66.1 ± 3. SNU 

64.1 SNU 

67.1 ± 4.4 SNU 



BOREXINO ^Be 

Standard model 
Seismic model 


3.36 ± 0.36 lO^cm 2s 1 
4.72 lO^cm-^s-i 3.045 ± 0.35 lO^cm-^g-i 


Water detectors 
SNO 

SNO +SK 

Standard model 2009 
Seismic model 


Predictions or Detections electronic neutrino flux 
5.045 ± 0.13 (stat)± 0.13 (syst) lO^cm-^s-^ 
5.27 ± 0.27 (stat)± 0.38 (syst) lO^cm 1 

4.21 ± 1.2 lO^cm-V^ 
5.31 ± 0.6 lO^cm-^s"^ 


neutrino flux electronic + other flavors in lO^cm ^s ^ 
SKI (5 MeV) 2.35 ± 0.02 (stat) ± 0.08 (syst) 
SNO D2O (5 MeV) 2.39 ± 0.23 (stat) ± 0.12 (syst) 
BOREXINO (2.8 MeV) 2.65 ± 0.44 (stat) ± 0.18 (syst) 



close to the detected neutrino fluxes, and it has been prominently used as a guide in 
the search for gravity modes (Turck-Chieze et al., 2004b, Garcia et al., 2007). It also 



Solar neutrinos, helio seismology and the solar internal dynamics 



50 



contributes to estimating the quality of predictions of the standard solar model. We 
consider that this seismic solar model is the only one for which we can reasonably put 
some error bars on the neutrino fluxes. These error bars come from the measurements of 
the different physical inputs introduced in stellar equations. Any error bar attributed to 
a standard solar model prediction can only be considered as a minimum error because it 
cannot be excluded that this model neglects important additional processes. See section 
6.5 for the oscillation parameters used. 



6. 5. Comparison with predictions including mixing parameters 



Clearly, the comparison of all these neutrino detection results with the predictions of the 
seismic solar model or the standard solar model demonstrates the existence of the MSW 
effect (section 3). Figure 9 shows the energy dependence of the neutrino oscillation 
due to matter or vacuum effects. Again, this is a very important result produced 
by the solar neutrino community. The predictions of Table 9 include the values of 
Bahcall fc Pena-Gar"ay| ( |2004[ ) of Am^2 = 710"^eV^ and tg'^612 = 0.45 (neutrino mass 
differences and mixing angles). All the predictions agree well with the detection results, 
and this agreement is much better than the one obtained from the standard solar model 



including the new updated CNO abundances (Turck-Chieze et al. 2004a Bahcall et al 



2006). Figure 10 (Aharmim et al. , 2008) summarizes how the oscillation parameters 



for the neutrino mass-eigenstate couple (1,2) are more and more constrained by the 
different results, and how these results lead (independently of solar models) to the 
values: Amj^ = 7.5910-5eV2 and t^^^is = 0.468. 



7. Beyond the standard solar model: a dynamical view of the Sun 



The equations used in the SSM framework (sections 4-5) do not consider the internal 
rotation and magnetic field. Twenty years ago, such a restriction was justified and 
did not prevent the description of the main phases of stellar evolution. However 
the description of the early and late stellar-evolution stages needs to be furthered. 



Young stars sustain a very strong magnetic field (Feigelson, 2003) and massive stars 



are strongly deformed by centrifugal forces and rapid rotation up to 50% like in 



the case of Achernar (de Souza, Kervella & Jankov 2003; Skumanich, MacGregor & 



Jackson, 2004). Moreover, the Sun is clearly a magnetic star (see the variability of the 



fundamental quantities given in Table 3) and its activity impacts the Earth environment. 
Because the origin of solar magnetism is internal, it is necessary to introduce the effects 
connected to the internal rotation and magnetic fields in solar models, to first replace 
the SSM, then the SeSM, by a more realistic DSM (Dynamical Solar Model). The most 
recent observations obtained with the SoHO satellite and the ground-based experiments 
contribute to the development of this new field. The Sun is rotating differentially at the 
surface, in roughly 24 days at the equator and 30 days at the poles. The internal rotation 
has a well known effect on the dynamo process which is characterized by variable activity 
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Figure 11. Top: Comparison between the solar internal rotation profile 
predicted by different solar models and the one deduced from helioseismology. 
From Turck-Chieze et al. (2010a). For models B5 and C^, the plotted value is 



ri/27r-0.2325 rather than r2/27r in order for the surface value they reach to match 
the one reached by model Be. The data points down to r/R© = 0.2 are deduced 
from the acoustic mode splittings determined by the observations of the GOLF, 



MDI and GONG instruments, from EfF-Darwich et al. (2008); Garcia et al. 



(2011). The data points in the core correspond to the core rotation extracted 



from the potential gravity modes observed. Bottom: Potential configuration of 
a fossil field in the radiative zone of the Sun: left Isocontours of the poloidal 
magnetic field Bp(r, 0) in G in meridional cut. c) Isocontours of the azimuthal 



magnetic field B0(r,0) in G in meridional cut. From Duez, Mathis & Turck- 



Chieze (2010) 



cycles of about 11 years with large fluctuations in intensity and duration. Therefore the 
knowledge of the solar internal rotation profile is a crucial ingredient of the dynamics 
of the Sun which will help to understand the solar activity. 
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7.1. The internal rotation and its consequences 



The Sun is the only star for which we can access the internal rotation in detail, thanks 
to the million of oscillation modes (million if we treat the different azimuthal orders m 
as separate) observed at its surface. The detection of the different components of the 
observed oscillation modes allows the determination of the solar interior rotation profile. 



1 

27 



^0 



Kn,^,^(r,^)fi(r,^)drd^ + en, 



7.1 



7.1.1. Rotation profile in the radiative zone Until recently, the only motion included in 
SSM was the gravitational settling (GS) (see section 4.3.4). This phenomenon describes 
the temporal migration of the heavy elements in comparison with the lighter ones toward 
the solar center. It is a very slow process, and equation (4.19) uses a diffusivity coefficient 
Di ~ 10 cm^/s. For comparison, the turbulent diffusivity Dt reaches 10^ cm^/s at the 



base of the convective zone (Brun, Turck-Chieze & Zahn 1999). GS and turbulence 



operate during the entire life of the Sun and consequently the current solar photospheric 
composition differs from its initial value by about 10%. Helioseismology, by measuring 
an helium content at 0.9 R0 of 0.243 that differs from the initial helium of the model 
of 0.277, validates this process. See table 6 for the impact of GS on the neutrino 
predictions. 

The rotation profile in the radiative zone has been difficult to establish. This is 
partly due to the small number of modes that reach this region (acoustic modes of 
degree C < 40) and the correction that needs to be applied to reduce the influence 
of the external differential rotation. Figure 11 (top panel) summarizes some current 
estimates. The detection of asymptotic dipolar gravity modes (Garcia et al. , 2007) 



favors, like the first gravity-mode candidates of Turck-Chieze et al. (2004b), an increase 



of the rotation rate in the core. The very central region seems to rotate faster than 
the rest of the radiative zone by between 5 to 8 times the corresponding solid rotation 



rate of 430 nHz (Eff-Darwich et al. , 2008). 430 nHz corresponds to the rotation rate of 



the convective zone at a latitude of 40 degrees. The transition between the radiative 
and convective zones generates an horizontal turbulent flow in a narrow region. Such 
hydrodynamic phenomenon occuring in a transition region from differential to solid- 



body rotation gave the name of tachocline to the corresponding layers (Spiegel & Zahn 



1992). We studied in great details (Turck-Chieze et al. 2010a) the transport of angular 



momentum induced by the rotation, by following the evolution of the Sun since the 
initial phase of contraction and introducing in the angular-momentum equation a term 



of advection and a term of diffusion (Zahn, 1992 Mathis & Zahn, 2005). We note that 



the rotation profile is mainly established during the early evolutionary phase and the 
present rotation rate favors the idea that the young Sun was a rather slow rotator. We 
also observe that there is a meridional flow in the radiative zone of the present Sun, but 
with a very low velocity of the order of 10~^ 10~^ cm/s, unlike the velocities observed at 
the top of the convective zone (order of a few m/s). This naturally produced interface 
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layers (tachocline) which are both very turbulent and dominated by horizontal velocity. 
It was also shown that the presence of such a rotation profile has almost no impact on 
the sound-speed profile, but certainly sustains an internal deep magnetic field. 



7.1.2. Meridional Circulation in the convective zone It is known since Galileo Galilei 
that the Sun rotates differentially in latitude and that the equator rotates faster by 30% 
than the poles, partly due to the thermal instability generated by the global rotation. 
Moreover it has been measured that such differential rotation is maintained in the 



entire convective zone (Kosovichev et al. 1997 Thompson et al. 2003). The meridional 



circulation (MC) of a large-scale flow velocity of 10-15 m/s, is observed at the solar 



surface and directed from the equator to the poles (Duvall, 1979). A return flow (still 
unobserved) must take place deeper in the convective zone, from the poles to the equator. 
This MC is probably playing an important role in the generation of magnetic field by the 



dynamo effect {e.g. (Brun & Rempel, 2009). Indeed, it transports not only matter but 



also magnetic flux (and angular momentum) with a timescale close to the solar cycle 



period (11 years), a crucial ingredient in fiux-transport dynamo models (Choudhuri, 
Schussler & Dikpati, 1995). Local helioseismology is conveniently used to invert the 



flow amplitude and direction in the upper solar layers, down to a few Mm below the 



surface. For instance, ring-diagram analysis (Haber et al. , 2002) reached a 15 Mm depth. 



Time-distance analysis (Zhao &; Kosovichev, 2004) confirmed their result but failed to 
detect any return flow. 



7.1.3. Zonal Flows in subsurface layers Zonal flows, or torsional oscillations (Howard 



& Labonte, 1980), are latitudinal bands of modulation of the rotational velocity. The 



mid-latitude bands move toward the equator during the activity cycle, while the high- 
latitude bands move toward the poles. The amplitude of these velocity modulations is of 
the order of 5 m/s. They are thought to be produced by feedback from the dynamo effect 



on differential rotation (Spruit, 2003). Kosovichev et al. (1997) detected the zonal flows 



in helioseismic data using surface-gravity waves or f-modes. Global helioseismology can 
only detect symmetrical patterns around the equator. It was shown that the torsional- 
oscillation pattern reaches deep inside the convective zone. Local helioseismology is 
needed to detect potential differences between the two solar hemispheres. It gives us 
access to these torsional oscillations and their depth dependence. The evolution of 



the zonal flow pattern during the solar minimum is being studied (Howe et al. , 2009) 



and is now used to forecast the strength of the upcoming solar cycle. Indeed, unlike 
the butterfly diagram that was traditionally used to observe the evolution of the solar 
activity cycle, the torsional oscillations are clearly visible even in absence of sunspots 
(which do not always emerge). The solar magnetic activity never stops. 
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7.2. The internal magnetic field 

The interaction between internal magnetic field and neutrinos would be very interesting 



to study (Couvidat, Turck-Chieze & Kosovichev, 2003 Rashba et al. 2007), but this 



phenomenon is probably small in the solar interior. 

We believe that a fossil field generated during the early phases of the formation 
of the Sun, when it was totally convective, may have survived up to the present age. 
It is now well established that if such a field still exists, it diffused over a timescale of 



Gyrs and must be constituted of a mixture of poloidal and toroidal fields (Zahn, Brun 



& Mathis, 2007 Duez, Mathis & Turck-Chieze, 2010). If it exists, we must determine 



its amplitude and its role on the variability of the solar activity, and on the potential 
magnetic interaction with neutrinos. 

The discovery of the so-called tachocline region emphasized the urgent need for 
the solar community to better understand the magneto-hydrodynamics of the radiative 
interior, and of the tachocline, and their coupling to the upper solar layers. It now seems 
well established that the tachocline plays a crucial role in the solar dynamo, since it is 
most likely the layer where the mean toroidal magnetic field, thought to be at the origin 
of the surface sunspots and butterfly diagram, is streched, amplified (by at least a factor 
of 100) and stored until it becomes magnetically buoyant ( Parker[ 1993; Rempel 2003). 
However little is known about the dynamics of the solar tachocline: is it turbulent or 
laminar, what types of circulation are present, what is the dynamical influence of the 



magnetic fields, why is it so thin (extending at most over 5% of the solar radius)? Spiegel 



& Zahn ( 1992 ) were the first to address directly some of those questions but in the purely 
hydrodynamical context. They showed that if no process were present to oppose its 
radiative spread, the solar tachocline would extend over 30% of the solar radius after 4.6 
Gyr, in complete contradiction with current helioseismic inversions. They demonstrated 
that the anisotropic turbulence could hinder the spread of the solar tachocline to only 



a few % of the solar radius. Elliott (1997) confirmed their results numerically with 



a 2D axisymmetric hydrodynamic code. Miesch (2003), using a thin layer version of 



the ASH code, showed that the coupling between randomly forced turbulence and an 
imposed shear flow gives rise to Reynolds stresses that transport angular momentum 



such as to reduce the shear. However several authors 


Rudiger & Kitchatinov, 


1997 


Gough & Mclntyre 1998 MacGregor & Charbonneau 




1999 Garaud 


2002 


Garaud & 



Garaud , 2008) proposed that the magnetic torques exerted by a weak internal fossil 
magnetic field (if it exists) could oppose the inward thermal hyperdiffusion (or viscous 
diffusion when thermal effects are neglected) of the solar tachocline. Such models favor 



a slow, rather laminar version of the tachocline. Gilman and collaborators (Dikpati 



Gaily & Gilman, 2004) developed a series of models that showed that the tachocline 
could become unstable through magnetic instability of toroidal structures embedded 
within it, resulting in a latitudinal angular momentum transport that suppresses the 
shear and limits its inward diffusion. This increasing number of solar tachocline models 
demonstrates how important it is to characterize the dynamical properties of the solar 
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radiative interior. 

With the advent of powerful supercomputers, non-hnear studies of the solar 
radiative interior and tachocline in full 3D MHD simulations became possible, but 
they lack the important and more detailed observational constraints that are required 
to progress in our understanding of the magnetohydrodynamics of the solar radiative 
interior. A space project like DynaMICCS (?), constraining the solar radiative 
dynamics, is clearly required to support and guide the 3D simulations. 



7.2.1. Time-distance helio seismology to study subsurface activity Distinct branches of 
local helioseismology have emerged: time-distance analysis (Duvall et al. , 1993), ring- 
diagram analysis (Hill, 1988), and acoustic holography (Lindsey & Braun, 1990) are 
the most widely used techniques. Local helioseismology is especially conducive to the 
study of the dynamics of the solar upper layers. For instance, it allows the analysis 
of meridional circulation, zonal flows, emergence of active regions, supergranular flows, 
and other phenomena closely related to the generation of magnetic fields by dynamo 
effect in the solar interior. The potential impact of these large-scale magnetic fields 
on the propagation of solar neutrinos (through, for instance, the resonant spin-flavor 
precession process) makes local helioseismology directly relevant to this review. For the 
sake of brevity, we will only introduce the time-distance analysis. 

Time-distance analysis measures the travel times of wavepackets propagating 
between different points on the solar surface using a cross-covariance technique. These 
wavepackets can be acoustic or surface-gravity waves. If the conditions on the 
photosphere were similar to the ones encountered on the Earth crust, it would be 
possible to follow the temporal evolution of a wave excited by a "sunquake" and measure 
its relevant travel time(s) without elaborate mathematical techniques. However, the 
perpetual stochastic excitation of waves by millions of convective cells of the upper solar 
layers makes any attempt at directly measuring any travel time of a specific wavepacket 
quite challenging. Hence the need to compute cross-covariances. 

We shall review a few results relevant to the potential interaction of neutrinos 
with the solar magnetic fields, because these results can (and are) used to discriminate 
between different theories of the formation of large-scale fields in the Sun. 



7.2.2. Structure and Dynamics of Sunspots Sunspots are regions of the solar surface 
that appear darker than their surrounding due to a lower temperature (see, e.g., the 
review by Solanki (2003). They are regions of intense magnetic fields (up to about 
3000 Gauss in the umbra) and are directly connected to the generation of a toroidal 
field by the dynamo effect. Although sunspots have been known for centuries, many 
details regarding their nature are still elusive. They are thought to be the consequence 
of a magnetic-field fiux raised by magnetic buyoancy (the density and pressure inside a 
fiux tube are lower than the surrounding plasma because the magnetic pressure adds to 
the gas pressure to satisfy the hydrostatic equilibrium), and piercing the solar surface. 
Local helioseismology proved a powerful tool to investigate their subsurface structure. 
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The first inversion of tlie sound-speed profile below an active region was presented in 



Kosovichev, Duvall & Scherrer (2000). Since then, numerous articles have dealt with 
this topic {e.g. Jensen et ai, 2001; Couvidat et ai, 2004). The main result is the 
existence of a two-region structure: a decrease in the sound speed compared to the 
quiet Sun, immediately below the solar surface, followed by an increase in sound speed in 
deeper layers. The reality of this two-region structure is controversial, due to limitations 
and approximations made during the inversion process. In particular, the impact of a 
magnetic field on the propagation and conversion (into fast MHD waves) of acoustic 



waves (Crouch & Cally, 2005) is neglected. Similarly Zhao, Kosovichev & Duvall (2001) 



applied the time-distance formalism to derive the flow velocities underneath a sunspot 
(down to 30 Mm), and found results consistent with the Parker model of sunspots 
(cluster model). Knowing the structure of sunspots can shed light on the interaction 
of the magnetic field fluxes with solar convection (strong magnetic fields are thought 
to impede convection), on the way these flux tubes arise and decay, where they are 
generated and anchored, and so on. 

7.2.3. Supergranulation as a Travelling Wave Using surface-gravity waves to measure 
the horizontal divergence of the flows in supergranules (convective cells of typical size 30 



Mm) Gizon, Duvall & Schou (2003) found that supergranulation undergoes oscillations 
and supports waves with periods of six to nine days. These waves are predominantly 
prograde, which explains the apparent super-rotation of supergranules (supergranules 
appear to rotate faster than magnetic features at the solar surface). The nature of 
this wave-like behaviour is not understood, but numerical simulations of convection 
in oblique magnetic fields showed the existence of traveling waves. Therefore, local 
helioseismology studies of supergranulation could shed light on the subsurface solar 



magnetic field. Indeed, Gizon, Duvall & Schou (2003) suggest that their result could be 



the beginning of the use of supergranular waves to probe the upper convection zone. 



7. 3. Toward a dynamical model of the Sun 

Building a dynamical model of the Sun is a difficult task, it is why we use today the 
seismic model to predict the different observed quantitites. The Dynamical Solar Model 
takes into account all the dynamical processes: the present magnetic variability of the 
external layers together with the slow evolution of the activity along the solar life. 
The history of the momentum transport will include the time evolution of a potential 
fossil field (introduced to describe the very young solar analogs) in interaction with 
the dynamo field and a good understanding of the phase of accretion ejection which 
may have modified the initial solar mass and its impact on the formation of planets 



(Turck-Chieze, Piau & Couvidat 2011). This model will include also the role of the 



gravity waves along time. It could take 10 years to build it with a lot of consequences 
on different fields of physics so, building such model justifies the instrumental progress 
one can perform during the next decades. 
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8. Conclusion, Open questions & Perspectives 



The solar neutrino puzzle stimulated impressive instrumental developments on both 
the physics and astrophysics sides. This puzzle was solved by proving the existence 
of the neutrino flavor oscillations, thanks mainly to the SNO detector but also to 
SuperKamiokande and BOREXINO. These results help us put strong constraints on 
the oscillation parameters 612 and 5mi2 and on the survival probability Pee(E). All 
the neutrino detectors, including the gallium ones, show the deflcit and support the 
LMA neutrino-oscillation solution. In parallel, the seismic investigation of the solar 



core with the very convincing GOLF and MDI instruments (Turck-Chieze et al. 2001b 



?), recently conflrmed by the ground-based networks GONG and BiSON (Basu et al. 



2009), shows the reliability of the neutrino flux predictions deduced from models taking 
into account the observed solar sound-speed proflle confronted to the actual neutrino 
detections, on the contrary the SSM predictions vary with time and appear now in 
marginal agreement with all the experiments. In this review we focused on the main 
facts demonstrating how successfully this complex problem was tackled. 

Of course, in light of this success, new dramatic questions have been raised and 
agreement does not mean complete understanding of the solar interior. It is not totally 
estabhshed if the SSM is the appropriate picture or if dynamical effects play some crucial 
roles in the context described in this review. We have shown the different directions 
of improvements which require to take into account all the components of the solar 
magnetism and the development of new laboratory facilities which produce equivalent 
plasmas to check their microscopic properties. 



8.1. Secondary effects of the neutrino properties 

Variability has been searched for in the neutrino detections because they may reveal 
other properties of the neutrinos: day-night variability can unveil some oscillation 
characteristics, and longer time variability can be the signature of some sensitivity to the 
magnetic moment of the neutrino. The global flt to the SNO day - night energy spectra, 
plus data from other solar neutrino experiments, strongly favored the LMA solution in a 



2-flavor MSW neutrino oscillation analysis (Ahmad et al. , 2002). Long term variability 



has not been convincingly evidenced in the SNO and SK detectors (Aharmim et al. 



2010a). It has been thoroughly investigated in the chlorine and gallium experiments 



(Gavryuseva &; Gavryusev, 


2001 


Sturrock & Scargle 


not been put deflnitively in evidence ( 


Sturrock 


2009 


)• 



time experiments, the statistics is high but the solar neutrino counts are the result of a 
complex analysis and of a strong reduction of the initial statistics; (2) in radiochemical 
experiments the statistics is, on the contrary, very low, and some differences in the 
extraction procedure over time may have perturbed the time analysis of the signal; (3) 
the magnetic moment is probably too small. The upper limit is presently less than 
10 ~^^fiB with BOREXINO and the internal solar magnetic fleld could be too weak. 



probably smaller than 3 MG even in the radiative zone. (Duez, Mathis & Turck-Chieze 
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Figure 12. Order of magnitude of the leptons and quarks. From Giunti & 



Studenikin (2009). 



2010) Therefore, a clear signature of magnetic interaction seems difficult to achieve. 



Nevertheless, more and more works are dedicated to this kind of interaction for the Sun 



with the hope of detection of a signal in BOREXINO (Das, Puhdo & Picariello, 2009), 



and with important consequences for more massive stars (Heger et al. , 2009). 



8.2. The neutrino masses 

Solar neutrinos contributed to develop neutrino physics. The sensitivity to sin^ 6*13 (la 
error) of the solar neutrino detectors (including BOREXINO, LENS, SNO+) can be 
improved down to 0.01 - 0.02 by precise measurements of the pp-neutrino flux and 
the CC/NC ratio, as well as spectrum distortion at high energies. The combination 
of experimental results sensitive to the low and high energy parts of the solar neutrino 



spectrum resolves the degeneracy of angles 6*13 and 6*12 (Goswami & Smirnov, 2005). The 



comparison of the present sin ^13 = 0.017 ± 0.026 as well as sin 612 measured in the 
solar neutrinos and in the reactor /accelerator experiments may reveal new effects which 
cannot be seen otherwise. This field has already benefited from atmospheric neutrinos 
and reactors like KAMLAND (see its role in Figure 10). For the global analysis of 



the experiments, see Fogli et al. (2008 2010), and Gonzalez- Garcia, Maltoni & Salvado 



(2011) who include also the cosmological results. 



Then come Double Chooz and the first long baseline with SK, ready to improve 
the present situation. Moreover the double /3 decay at Mainz and Troitsk, MiniBooNE 
and MINOS add very important constraints on the effective electron neutrino mass 
which is now smaller than 1 eV. A lot of discussions appears around the potential 



existence of a fourth neutrino, the sterile neutrinos (Mention et al. , 2011) and other 



short baseline electron neutrino disappearance (Giunti & Laveder, 2009). If the neutrino 



were massless, neutrino and antineutrino would unambiguously be the same particle. 
However neutrinos have a mass, which results in two possibilities: neutrinos are a Dirac- 
type particle (neutrinos and antiparticles differ by a leptonic internal charge); or they are 
Majorana particles, and the neutrinos are the self conjugate of the antineutrinos. Today 
we know for sure that neutrinos are massive, see Figure 12 to see the order of magnitude 



of their fundamental constituents. See also Petcov (2009); Giunti & Studenikin (2009); 



Valle (2010) for general reviews on the present status of the field. The Large Hadron 



Collider, which just started running, could contribute to this success story by shedding 
light on the Majorana character of the neutrino. 
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8.3. Future experimentations: Neutrino Astronomy from Sun to Supernovae 

The currently running solar neutrino experiments and the upcoming ones (ICARUS, 



Arneodo et al. (2006), MiniLENS - LENS Raghavan et al. (2008)) are extending the 



energy threshold of neutrino detection toward the low energy (more statistics and 
less dependence on the theoretical spectrum for the boron neutrinos). Therefore 
they will be able to separate the pp contribution from the other neutrino sources, 
because the boron and beryllium neutrinos are measured independently; they will also 
allow a detection of the neutrinos coming from the CNO cycles, thus heralding the 
beginning of neutrino astronomy. The liquid scintillator of the SNO+ experiment in 
the SNOLAB (located in the same Sudbury mine than SNO, but at a 2 km depth), 
containing dissolved neodymium, will measure neutrinoless double beta decay in that 
isotope to check whether or not the neutrino is a Majorana particle. All the current 
underground laboratories will also be able to detect supernova explosions. Moreover, 
a network of neutrino detectors is now being installed around the world. If a burst 
of neutrinos is detected, SNEWS (Supernova Early Warning System) will alert the 
astronomy community (and the amateur astronomer community as well). 

On the astronomical side, the activity of the Sun is now at the center of the research 
effort, with the inclusion of all the different dynamical processes in solar models in 1, 2, 
or 3D. A new generation of missions are beginning to operate this year, SDO has been 



launched in February 2010 (Scherrer et al., 2002), PICARD in June 2010 (Thuillier, 



Dewitte fc Schmutz 2007). A new instrument GOLF-NG (Turck-Chieze et al., 2006 



2008 , 2010b ) has been qualified on ground to pursuit in space the gravity mode detection 



by a multichannel detection and add also some solar atmospheric investigation in an 
extended mission of formation flying (Turck-Chieze et al. 2009a) . The main goal of 



these missions is to better understand the Sun-Earth relationship in near real time (on 
a scale of days or less for the forecast of the Space Weather) or on a longer timescale 
(decade, to understand the space climate). Of course, this detailed information will be 



accompanied by the strong development of asteroseismology with the COROT (Baglin 



et al., 2006) and KEPLER (Koch, Borucki et al. 2010; Chaplin et al. 2010) missions. 



which will observe thousands to millions of pulsating stars. This dramatic evolution 
will give a new orientation to the Neutrino-Astronomic community: the detailed study 
of the explosion of different types of stars is bound to produce, once more, excitating 
times. 



8.4. Solar core and Dark Matter 

Despite the great success of the solar physics story, the central solar core below 0.10 R0 
containing about a quarter of the solar mass, is still not completely understood. Most of 
the acoustic modes have been detected but they do not allow us to properly describe the 
thermodynamics of this region of the Sun through the sound speed. Consequently the 
central temperature deduced from the seismic model is obtained by assuming that the 
temperature and density profiles follow the classical equations of stellar evolution. This 
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hypothesis leads to a good agreement between prediction and detected boron neutrino 



fluxes but it has been known for 20 years ( Giraud-Heraud et aL 1990 Dearborn et al. 



1990 Kaplan et al. , 1991 ) that these profiles could be modified by the presence of dark 



matter. This matter exists and must be present in the stellar cores. It has not yet 
been identified but it shall change the transport of energy (Spergel& Press, 1985) and 
this transport depends on the scattering cross-sections between dark matter and stellar 
plasma ( |Giraud-Heraud et al.[ |1990[ [Lopes, Bertone fc Silkf |2002t [Lopes fc Silk| [2002[ ). 
This matter thermalizes the central core and, due to the fact that the temperature stays 
relatively constant in this region, it modifies the asymptotic behaviour of the gravity 



mode frequencies. This could be an interesting signature to look for it (Lopes & Silk 



2010). The first detection of dipole gravity modes with GOLF/SoHO excludes a large 



range of WIMPs with mass < 10 GeV (Turck-Chieze et al. 2011b). This interesting 
investigation is an additional argument between others (Turck-Chieze et al. , 2005 , 2009a) 
to pursue the gravity mode detections with the new generation of instruments mentioned 
above. 



8.5. Last remarks 

The present ongoing experiments and others in construction shows the real beginning 
of Neutrino Astronomy. The detection of CNO neutrinos and pure pp will add new 
constraints on the knowledge of the central CNO in addition to the present revised 
photospheric values. The on going opacity measurements will add crucial checks to 
separate the effects of microscopic physics from macroscopic physics in the radiative 
zone. The g-modes will reveal the order of magnitude of the fossil field and its 
interconnection with the dynamo field, this progress will open a new page on the 
knowledge of the variability of our star. Double Chooz and Minos will constrain ^13 
which could reveal some difference between matter and vacuum. 

These last three decades have been extremely productive. It was not possible in this 
review to give credit to all the works devoted to this subject. This review is certainly 
not totally exhaustive but we describe the main highlights and the direction of new 
works. We limited the number of tables, so we encourage the reader to return to our 
first review of 1993 and to the tables published by J. Bahcall and his collaborators. 
We hope, for the youngest members of our community, that stars (including Sun) and 
supernovae will trigger studies as exciting as the ones summarized in this review. 
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